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ANALYTICAL INVESTIGATION OF SUPERSONIC TURBOMACHINERY BLADING 
I - COMPUTER PROGRAM FOR BLADING DESIGN 
by Louis J. Goldman and Vincent J. Scullin 
Lewis Research Center 

SUMMARY 

A FORTRAN IV computer program for the design of supersonic blading based on es- 
tablishing vortex flow within the blade passage is presented. The method of character- 
istics, as applied to the two-dimensional isentropic flow of a perfect gas, was utilized 
for the blade design. The equations necessary for the design are developed. The infor- 
mation required for the program consists of an inlet flow angle, specification of the in- 
let, outlet, and lower- and upper- circular-surface Mach numbers, and the specific-heat 
ratio. The program output consists of the blade coordinates and, if desired, a printer 
plot of the blade profile and flow passage. In addition, supersonic starting and flow sep- 
aration calculations are performed by the program and obtained as output. An example 
is included to indicate the use of the program and the results obtainable. 


INTRODUCTION 

Supersonic compressors and turbines are employed in special circumstances be- 
cause of their simplicity and low weight. A recent application for a supersonic turbine 
involves the hydrogen-fueled open- cycle auxiliary space power system descriDed in ref- 
erence 1. If the highest practical efficiency is to be obtained from supersonic compres- 
sors or turbines, proper design methods must be available. 

A method for designing supersonic blade sections based on two-dimensional isen- 
tropic flow is given in reference 2. The method consists of converting the uniform par- 
allel flow at the blade inlet into a vortex flow field, turning the vortex flow, and recon- 
verting to a uniform parallel flow at the blade exit. The application of this design pro- 
cedure involves specification of the inlet and outlet Mach numbers, the lower- (or con- 
cave) surface Mach number, the upper- (or convex) surface Mach number, the inlet flow 



angle, and the specific-heat ratio of the working fluid. In general, a wide range of de- 
signs is possible by selection of these parameters. Guidance in the selection of a blade 
design is obtained by considering blade shape, solidity, and supersonic starting and flow 
separation problems. In reference 2 the effect of some of the design parameters for low 
Mach numbers and a specific-heat ratio of 1. 4 is examined. 

In view of the interest in hydrogen-fueled auxiliary space power systems, an analy- 
sis was conducted to gain a better understanding of the effects of the design parameters 
on the resulting blade geometry and to extend the results of reference 2 to levels of in- 
terest for such systems. In reference 3, the effect of surface Mach numbers, inlet flow 
angle, and specific-heat ratio on the geometric characteristics of supersonic impulse 
turbine- blade sections is investigated over an inlet Mach number range of 1. 5 to 5. 0. 
Blade design limitations resulting from supersonic starting and flow separation problems 
are also considered. In the present report, a description and a FORTRAN IV listing of 
a computer program for the design of blading applicable for any supersonic Mach number 
level and specific-heat ratio are presented. Supersonic starting and flow separation cal- 
culations are also performed by the program. An example is included to indicate the 
use of the program and the results obtainable. The report is organized so that those 
persons desiring to use the program need only read the sections METHOD OF ANALY- 
SIS, DESCRIPTION OF INPUT, and DESCRIPTION OF OUTPUT. All necessary infor- 
mation pertaining to the program itself is contained in the sections DESCRIPTION OF 
INPUT, DESCRIPTION OF OUTPUT, and PROGRAM DESCRIPTION. 

SYMBOLS 


A area, ft 2 (m 2 ) 

a speed of sound, ft/sec (m/sec) 

C reduction in maximum weight flow due to two-dimensional flow (eq. (34b)) 

C dimensionless blade chord, chord/r* 
f(R*) function defined by eq. (10b) 

G dimensionless blade spacing, spacing/r* 

g c dimensional conversion constant, 32. 17 ft-lb/(lb)(sec 2 ) (1 kg-m/(N)(sec 2 )) 

h blade height, ft (m) 

j index for upper surface of blade 

* 

K dimensionless vortex constant defined by eq. (23) 
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K* value of K* for which weight flow is maximum 

max 

k index for lower surface of blade 

M Mach number, V/a 

M* demensionless velocity or critical velocity ratio, V/V cr 

(M*) maximum inlet velocity ratio for supersonic starting 

1 max 

(M*) minimum lower-surface velocity ratio from separation criterion 

^ min 

(M*) maximum upper-surface velocity ratio from separation criterion 
u max 


m slope of Mach line 

m slope of wall segment 

2 2 

p static pressure, lb/ft (N/m ) 

(p ) maximum lower-surface static pressure from separation criterion 

L max 

Q vortex flow parameter (eq. 34a)) 

R radius in vortex field, ft (m) 

R* dimensionless radius in vortex field, R/r* 

r* radius of sonic velocity streamline in vortex field, ft (m) 

V velocity, ft/sec (m/sec) 

V cr critical velocity, ft/sec (m/sec) 

w weight flow, lb/sec (kg/sec) 

w maximum weight flow, lb/sec (kg/sec) 

max 

x* dimensionless X- coordinate of blade (fig. 1), X/r* 

x* dimensionless x- coordinate of transition arc (fig. 1), x/r* 

Y* dimensionless Y- coordinate of blade (fig. 1), Y/r* 

y* dimensionless y- coordinate of transition arc (fig. 1), y/r* 

a circular arc turning angle, rad 

/3 flow angle outside blade passage, rad 

y specific- heat ratio 

0 total flow turning angle, rad 

p Mach angle, rad 



v Prandtl-Meyer angle, angle through which flow must turn from Mach 1 to re- 

quired Mach number, rad 

incremental flow turning, rad 

(";) maximum inlet Prandtl-Meyer angle for supersonic starting rad 

minimum lower-surface Prandtl-Meyer angle from separation criterion, rad 
(^ u ) maximum upper-surface Prandtl-Meyer angle from separation criterion rad 

ITlcLX 

P density, lb/ft 3 (kg/m 3 ) 

a blade solidity, C*/G* 

<P velocity direction angle, rad 

Subscripts: 

d downstream of normal shock 

i blade inlet 

j index for upper surface of blade 

k index for lower surface of blade 

l lower surface of blade 

max maximum 

min minimum 

o blade outlet 

u upper surface of blade 

Superscripts: 

total- state conditions 

METHOD OF ANALYSIS 
Blade Description 

The design of supersonic blade sections described herein is based on establishing 
vortex flow within the blade passage by a procedure analogous to that given in refer- 
ence 2. The blade so designed consists essentially of three major parts: (1) inlet tran- 
sition arcs, (2) circular arcs, and (3) outlet transition arcs. A typical blade, with per- 
tinent nomenclature noted, is shown schematically in figure 1. The inlet transition arcs 
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AB Inlet upper transition arc 

FG Outlet upper transition arc 

CD Inlet lower transition arc 

HI Outlet lower transition arc 


AF Upper circular arc Y* 



Figure 1. - Typical supersonic blade section. (All coordinates are made dimensionless by 
dividing by r*. ) 

(lower and upper) are required to convert the assumed uniform parallel flow at the blade 
inlet into vortex flow. The concentric circular arcs turn and maintain the vortex flow. 
Finally, the outlet transition arcs reconvert the vortex flow into uniform parallel flow at 
the blade exit. Straight-line segments parallel to the inlet and outlet flow directions 
complete the blade profile. Methods for creating a finite thickness at the leading and 
trailing edges are given in reference 2. 

In general, the inlet lower transition arc reduces the Mach number from its value 
at the blade inlet M^ to a preselected value of the lower-surface Mach number M^, 
whereas the inlet upper transition arc increases the Mach number to a preselected value 
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Figure 2. - Surface Mach number variation for typical 
blade section. 


of the upper- surface Mach number M^. The surface Mach numbers remain constant, at 
these preselected values, on the lower and upper circular arcs. At the outlet region the 
procedure is reversed. The surface Mach number variation is shown in figure 2 for a 
typical blade. 

The amount of flow turning produced by either the lower or upper surface of the 
blade consists, in general, of two parts (fig. 1): (1) the turning produced by the transi- 
tion arcs and (2) the turning produced by the circular arcs. When isentropic flow turn- 
ing at supersonic speeds is considered, it is convenient to introduce the Prandtl-Meyer 
angle v, which is defined as the angle through which the flow must turn from Mach 1 to 
the required Mach number. The flow turning produced by a transition arc is then equal 
to differences in Prandtl-Meyer angles and is ja - i> l and v Q - v, for the inlet and 
outlet lower transition arcs and and - v Q for the inlet and outlet upper tran- 

sition arcs, respectively. The turning produced by the inlet or outlet transition arcs 
cannot exceed the inlet or outlet flow angle ^ or /3 respectively. The relation be- 
tween Prandtl-Meyer angle v and Mach number M is given by the following equations 
(ref. 4): 
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( 1 ) 


!/=I 


4 




arc sin 




where the dimensionless velocity or critical velocity ratio M* is 


M* = 



( 2 ) 


This relation is shown in figure 3 for specific-heat ratios of 1. 3, 1. 4, and 1. 66. As the 
Mach number approaches infinity, v approaches an upper limit of | + l)/(y-l)-l]. 
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Blade Design 


The method of characteristics as applied to the two-dimensional isentropic flow of a 
perfect gas is utilized in the design of the supersonic blade sections. A description of 
the method of characteristics is given in references 4 and 5, and its application to the 
design of supersonic blade sections is given in references 2 and 6. For purposes of cal- 
culation, the flow field is considered to be divided into small regions, in each of which 
the flow properties are assumed to be constant. If adjacent regions are to differ slightly 
in properties, then the boundary between the regions must be characteristic lines and 
can also be shown (ref. 5) to be Mach lines. Therefore, each region is, in general, 
bounded either by a Mach line or a physical boundary. In figure 4 the flow field for a 


AB Upper transition arc 

CD Lower transition arc 

CE Major expansion characteristic 

AE Major compression characteristic 



Figure 4. - Characteristic network within blade passage. 


typical blade passage is divided by characteristic lines into a finite number of regions. 
The vortex-flow region is bounded by the circular arcs and the outermost vortex charac- 
teristics AE and CE. The transition arcs are composed of straight-line segments, and 
within each region bounded by these segments the flow is constrained to follow the wall 
direction. The mathematical equations necessary to define the blade are developed in 
the following sections. 

Circular arcs . - Within the concentric circular arcs, vortex flow exists; therefore, 
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VR = Constant 


(3) 


where V is the velocity and R is the radius in the vortex field. In this report, dimen- 
sionless parameters are used whenever possible; if this procedure is followed, equa- 
tion (3) can be rewritten as 



Constant 

V r* 
cr 


(4) 


where V is the critical velocity and r* is the radius of the sonic velocity stream- 
line in thevortex field. At R = r*, V = V cr ; therefore, the constant is V^r*. Equa- 
tion (4) then becomes 


M*R* = 1 


(5) 


where M* = V/V cr is the dimensionless velocity and R* = R/r* is the dimensionless 
radius in the vortex field. The Prandtl-Meyer angle v is related to M thiough equa- 
tion (1). Therefore, once v ^ and are specified, M* and M* are fixed, and the 
circular arc radii R* and R* are determined from equation (5). The amount of lower 
circular arc turning for the inlet and outlet portions of the olade ^ ^ and ot^ , Q , 
spectively, are 


(6a) 

and 

a l,o‘ fl o + K-''l ) (6b) 

Angles measured in the counterclockwise direction are considered positive. With this 
convention, the inlet flow angle 0. is positive, and the outlet flow angle 0 Q is negative. 
Similarly, for the upper circular arc 




(7a) 


and 


V 0 - + K - "o> 


(7b) 


The circular arcs are completely described by specification of v-, v Q , v v u , and 0^. 
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The outlet flow angle does not have to be specified because it can be related to M., 
M q , and /3. from the following consideration. 1 

Because the inlet and outlet blade spacing is the same, the inlet and outlet blade pas- 
sage areas A^ and A q , respectively, are related by geometry (see fig. 1) according to 
the equation 


A^ cos /L 

A o cos ^o 


The area ratio A^/A q can be obtained from the continuity equation (ref. 5) with the re 
suit that equation (8) becomes 


Cr 


P 0 = -arc cos-( 


M, 


M, 






cos /3 i ^ 


(9) 


Lower transition arcs . - The lower surface is composed of an inlet and outlet tran- 
sition arc and a circular arc (fig. 1). For symmetric blades (i. e. , t/. = v ), the two 
transition arcs are identical, and, therefore, only one needs to be calculated. For 
asymmetric blades, the two arcs are not identical, one being smaller (less turning) than 
the other. However, the smaller transition arc corresponds to a portion of the larger 
arc, and, consequently, only the larger arc need be calculated. If v. is greater than 
V transition arc is the larger of the two arcs. For simplicity, the inlet tran- 

sition arc is assumed to be the larger arc in the following discussion. 

In figure 5, the lower transition arc is shown, and the nomenclature used in the 
computer program is indicated. The calculations are performed with respect to the 
nondimensional axes x* and y* (where the x- and y- coordinates are made dimension- 
less by dividing by r*). The transition arc coordinates are generated in a sequential 
manner (starting at x* = 0, y* = R*) by obtaining the intersection of the straight-line 
wall segments and straight Mach lines for a specified small change in flow turning. The 
Mach lines are detei mined from the outermost or major vortex- expansion characteris- 
tic, and the wall segments are determined from the flow direction. After the transition 
arc coordinates are calculated, they are rotated through an angle of a ^ . to obtain the 
coordinates of interest in the blade design (see fig. 1). 

For the vortex region, it can be shown (ref. 2) that the velocity direction y and the 
dimensionless radius R* are related along the characteristic line by the equation 
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transition arc 



Figure 5. - Nomenclature used for calculation of inlet lower transition arc 


(10a) 


(p = ± — £(R*) + Constant 
2 


where 



(10b) 


Two families of characteristics exist: the positive sign in equation (10a) gives the ex 
pansion lines, and the negative sign gives the compression lines. The major vortex- 
expansion- characteristic equation is 


<p =|[f(R*) - f(R*)] 


( 11 ) 


since the boundary condition at x* = 0 is that cp = 0 and R* = R*. if the flow field is 
divided into small regions, the flow direction at any point along the characteristic line 
(as given by eq. (11)) may be considered to be equal to the flow direction within the adja- 
cent flow region, as indicated in figure 5. For k transition arc segments, each of 

which produces At/ degrees of turning, the flow direction within any flow segment w 
is given by k, i 

^k, i = v i " v i ' ( k " 1 ) Al> = <?k+l, i + A " (12) 

where k is an integer that varies from 1 to [(i/. - v^/Av] +1. At kJ(i/.-i/ )/Av 1 + 1 
the flow direction is 0, and at k = 1, it is ^ V Equating equations (11) and (12)' and ’ 
eliminating f(R*) through equations (1), (5), and (10b) result in 

f ( R k, i) = 2v \ - | " 2 ^ k " (13) 

which relates the dimensionless radius R* to the incremental flow turning Au along 
the major vortex-expansion characteristic. At k = [(y. - i^)/ai>] +1, R* . = r*. as k 
is decreased, R^ ^ decreases and is obtained by the simultaneous solution of equations 
(10b) and (13) by kn iterative procedure. Since R* . can be determined for any value of 
k, the x*, y* coordinates along the major expansion characteristic can be obtained from 
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x k, i = -K, i sin v k _ j 


(14a) 


and 



cos 


«v 


(14b) 


These coordinates are also points on the straight Mach lines. The equation specifying 
the Mach line is therefore determined once the slope is obtained, which is easily accom- 
plished because the Mach line is inclined at the Mach angle j u to the velocity direction 
(see fig. 5). The Mach line, therefore, makes an angle of <p + fi with respect to the 
x* axis. It is possible to define the Mach line at the mean Mach angle to the mean flow 
direction (ref. 5) so that the slope of the Mach line m^ . is given by 


m 


k, i 


tan 


^k,i 


^k+l,i } 4 k,i + M k+1,0 


(15a) 


where 


u, . - -arc sin 
k ’ x \M 


-arc sin 


k, i/ 


VP 


+ l\o*2 




' y 


(15b) 


The equation of the Mach line is therefore 


y* = m k, i< x * - x £, i) + y k, i 


(16) 


where k varies from 1 to k = (v- - v,)/kv. The equation of the transition arc seg- 

max 1 l 

ment is more easily determined since each segment is a straight line parallel to the ve- 
locity direction cp. The slope of the wall segment m^ . is then 


m k,i = tan ‘Pk+l.i 


(17) 


and the equation for the wall segment is 


y* = m k i X * " ( X 7*> + 

k ’ 1 1_ 1 k+1, ij L k+1, i 


( 18 ) 
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where k varies from 1 to k max and x* and y* are the lower transition arc coordi- 
nates. The values of x* and y* are known at k = ( v ■ - v,)/ Lv + 1 = k +1 where 
^ ^ ^ it max ? 

X* = 0 and y* = . The remaining transition coordinates are generated by finding the 

intersection of the Mach lines with the wall segments starting at k = k and sequen- 
tially decreasing k until k = 1. The intersection of the two straight lines is given by 


( x f) 

1 M 


(yp 


k+1, i 


- m 




k+1, i_ 




m, .x,* .) 
k, i k, v 


m 


k, i " m k, i 


(19a) 


and 


(yf) 

1 k,i 




(yf) 


k+1, i 


- m, 


k. i< x ?) 


- m k i(yjt ; 
k+1, ij k ’ 1 k » 1 


- 111 . 


k, i x k, P 


m k, i " m k, i 


(19b) 


The transition arc coordinates obtained from equation (19) are rotated through an angle 
a l,\ resultin g in the X*,Y* coordinates of interest in the blade design. The rotated 
coordinates X* and Y* are obtained from 

(X*) = (x*) cos a . - (y*) sin a . (20a) 

L k, i L k, i 1 > 1 4,i 1 > 1 


and 


(Y*) = (x*) sin a . + (y*) cos a, . (20b) 

1 k,i 1 k, i L > 1 4,i 1 > 1 

Upper transition arcs . - The upper surface (like the lower surface) is composed, in 
part, of an inlet and outlet transition arc, only one of which must be calculated. For 
simplicity, it is again assumed that the inlet transition arc is the larger of the two arcs. 
(For the upper arc, this requires that v Q be greater than v^) The upper transition 
arc is shown schematically in figure 6, and the pertinent nomenclature is noted. The 
procedure employed to calculate the upper transition arc is analogous to that used for 
the lower transition arc; the resulting equations which are of similar form are not re- 
peated herein. The subscript j is used to represent the upper transition arc coordi- 
nates where j varies from 1 to i = (u - v.)/Ai> 

Geometric parameters . - After the blade calculations have been performed, a num- 
ber of blade parameters of interest, including blade solidity, spacing, chord, and total 
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Major vortex-compression characteristic 



Figure 6. - Nomenclature used for calculation of inlet upper transition arc. 



flow turning angle, are calculated. The blade spacing G* and chord C* are obtained 
from the blade coordinates, and the solidity is obtained from the ratio C*/G*. The total 
flow turning angle 0 is obtained from the inlet and outlet flow angles. The blade co- 
01 dinates are also translated by G* so that coordinates for two complete blades are 
obtained. 


Design Limitations 

The design limitations (i. e. , the constraints on the choice of 1 / and v for spec- 
ified i/., v Q , /3p and y) imposed by consideration of supersonic starting and flow sepa- 
ration problems have been discussed in reference 3. These limitations are calculated 
by the procedure described in the following paragraphs and are given as output from the 
computer program presented herein. 

Supersonic starting . - The problem of establishing supersonic flow on startup is 
discussed in reference 7 for supersonic compressors. For supersonic turbines, the re- 
sulting design limitations due to starting are presented in reference 2, where it is as- 
sumed that a normal shock wave spans the blade inlet at the instant of startup. Under 
this condition it is necessary to ensure that the weight flow can pass through the turbine. 

The maximum value of the inlet Prandtl-Meyer angle (v.) is determined by first 

1 max 

finding the maximum weight flow through the blade passage, while taking into account the 
normal shock losses. This maximum weight flow is then equated to the flow rate after 
the shock has passed through the passage. 

The weight flow through the passage is obtained by integrating the continuity equation 
in the vortex region 


w = 



pV dR 


( 21 ) 


where p is the density, V is the velocity, and h is the blade height. The density can 
be written as (ref. 2) 



( 22 ) 


where ^ and a T .^ ^ are the density and the speed of sound, respectively, just down- 
stream of the shock, and are evaluated at total conditions. For perfect gases, the total 
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temperature is constant through a normal shock so that a'^ ^ - a 'p Utilizing equation (3) 
and the definition 



results in equation (21) in the following form: 



(24) 


Differentiating equation (24) with respect to K* and setting the result equal to 0 give the 
value of K* (denoted as K* iax ) for which the weight flow is a maximum. This proce- 
dure gives 



Changing the variable from R to M* gives, for equation (25), 
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~1 1/ (r- 1) 


M* 2 


dM* 

M* 



(26) 


Integrating the right side of equation (26) results in 



1-K 



^/(y-l) 


*2 

~ maxl m; 


(27) 


Similarly, equation (24) in terms of M* becomes 



The value of Kj^ is determined from equation (27) by an iterative procedure (for 
given values of M* and M*), and then the maximum weight flow is obtained from equa- 
tion (28). After the shock has passed through the turbine, the weight flow can again be 
obtained from equation (2 1), which is rewritten by substituting VR = V r* and chang- 
ing the variable from R to M* to give U 


w = 



(29) 
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Substituting the following relations (ref. 5) into equation (29) 




p = p 


1 \y + !/ 


( 30 ) 


and 


V = a 

v cr 




(31) 


gives 


w = r*ha. 


W^t-TT r” /,V 


i M* 2 V /(y ~ 1} dM * 


(32) 


M 


where the term multiplying r^hajp! is defined as the weight-flow parameter. Equating 
the two weight flows, equations (28) and (32), results in 


Q _ P i, d _ p i, d 


1 - C p 


(33) 


where 


»_ M * M u /’ M “ (z±l - nzJ. M * 2 ) 1/,<r l> — 

l 2 2 ) M* 


(34a) 


and 
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The parameter C has been shown in reference 2 to be the reduction in maximum flow 
rate caused by two-dimensional flow. The quantity y P ! is the total pressure recov 
ery for a normal shock and is given by (ref. 5) ’ 




2 y/(y- 1) 
max 



(3 5) 


The maximum inlet velocity ratio (M*) is obtained (lor given values of M* and 

M u ) by simultaneous solution of equations (33) and (35), by using the definition in equa- 
tions (34a) and (34b), and by obtaining K^ ax from equation (27). The maximum inlet 
Prandtl-Meyer angle for supersonic starting ( v ) is then obtained from equation (1) 

-CM A 

— w separation . - Analysis of supersonic blade sections (ref. 3) has shown the de- 
sirability of maintaining high surface Mach numbers to alleviate the problem of super- 
some starting. Under these conditions, however, adverse pressure gradients created 
on the blade surfaces would be expected to cause flow separation and, consequently, 
poor performance. Experimental investigation of simple shapes with incompressible 
flow at fairly high pressure gradients (ref. 8) has indicated that if the coefficient of 
pressure recovery (defined as the ratio of the pressure rise to the dynamic pressure at 
the initial point) is less than about 1/2, flow separation may be avoided. This criterion 
was used in the analysis presented in reference 3 to give some indication of the design 
restrictions due to flow separation. For supersonic velocities the separation value of 

the coefficient of pressure recovery may be less than 1/2. The calculation^ procedure 
is as follows. 

Flow separation can occur on both the lower and upper surfaces of the blade, but 
since the calculational procedure is similar for both cases only the derivation for the 
lower surface is presented. The flow separation criterion can be written as 


(P z ) 


max 

PiVf 

2g„ 


- Pi 


1 

2 


(36) 


where ( P/ )^^ is the maximum lower-surface pressure possible (for given inlet condi- 
tions) without causing separation. Equation (36) can be rewritten in the form 
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(P;) 

max 





(37) 


where the equation has been divided by the inlet total pressure p'j. Substituting the fol- 
lowing relations (ref. 9) into equation (37) 

-R = /l .ZJlIm* 2 Y /<V ’ 1) (38a) 

P V v + 1 / 

and 

J_ PV 2 _x_ m *2 

g c = r + 1 (38b) 

p i _ r :1 m* 2 

y + 1 


and simplifying result in 


r 


( M ?) 


min 




1 - 1-31 


y 


+ l 


! M* 2 


1 


JL-ml 


*2 


l + i 


+ 1 


, 1 - Z-li M* 2 


y + 1 


(y-l)/y 




1/2 


(39) 


Similarly, applying the same criterion to the upper surfaces gives 

/ - M M u ) 

1 / y + 1 l 




i - 




max 


1 


a -^ K ) 

y + i 



|(y-l)/y 


1/2 


> (40) 


The corresponding Prandtl-Meyer angles (v } ) and (v ) are obtained from equa- 

L min max 

tion (1). 


min 
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TABLE I. - INPUT FORMAT FOR FORTRAN COMPUTER PROGRAM 














DESCRIPTION OF INPUT 


The input for the computer program consists of an inlet flow angle, several Prandtl- 
Meyer angles, the specific-heat ratio y, and an angular increment. In addition, three 
optional switches must be set which regulate the content and form of the output. All the 
input parameters, except y and the switches, must be specified in degrees. The input 
variables are as follows: 

BETAN inlet flow angle, (3^ 

DELV flow turning increment (recommended value, 0. 1), Ay 
GAM specific- heat ratio, y 

IPRINT a value of 0 will result in printing of rotated blade coordinates X* and Y*; 

a value of 1 will cause both rotated and unrotated coordinates X*, Y* and 
x*, y* (see fig. 1) to be printed 

ISTART a value of 0 will cause both starting and blade design calculations to be printed 
out; a value of 1 will cause only starting calculations to be performed and 
printed out 

NPLOT a value of 0 will cause blade profile and flow passage to be plotted; a value of 
1 will suppress the plot 

VIN inlet Prandtl-Meyer angle, ja 

VLOW lower-surface Prandtl-Meyer angle, y ^ 

VOUT outlet Prandtl-Meyer angle, v Q 

VUP upper- surface Prandtl-Meyer angle, y y 

The flow turning increment Ay must be specified so that (y^ - y^)/ Ay, (y Q - y^)/Ay, 

- y.)/Ay, and (y y - y Q )/Ay are all integers. Table I shows a sample input card. 

DESCRIPTION OF OUTPUT 

An example of the output obtained from the program is shown in table II. The out- 
put corresponds to the input data shown in table I and consists of tables of coordinates 
for the description of two blade sections, supersonic starting and flow separation param- 
eters, and a plot of the blade profile and flow passage. This example required approxi- 
mately 0. 2 minute of computer running time. Each section of the output has been num- 
bered to correspond to the following description: 

(1) The first output of the program is a listing of the supersonic starting parameters. 
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If ISTART=0, the program will continue with the blade design calculations. If ISTART=1, 
no further output is ootained. 

(2) The next output is a listing of all the input data plus the value of the calculated 
outlet flow angle. 

(3) If IPRINT=1, the following tables are printed: 

(a) Uni otated coordinates of the lower- and upper-surface transition arcs 

(b) Coordinates of the lower- and upper-surface circular arcs 

(c) Coordinates of the upper- surface straight-line segments 

This output is, in general, not of interest except for debugging purposes and may be 
omitted by setting IPRINT=0. 

(4) The next output is tables of the rotated coordinates of the lower and upper tran- 
sition arcs. In addition, these coordinates are translated by the value of the blade 
spacing so that coordinates for two blade sections are obtained. Every tenth calculation 
point is printed if Av < 0. 2°; otherwise every calculation point is printed. 

(5) The next output is a listing of miscellaneous parameters including 

(a) Inlet, outlet, and surface dimensionless velocities or critical velocity ratios 

and Mach numbers 

(b) Dimensionless blade spacing, chord and solidity 

(c) Separation limitation for the lower- and upper-surface Prandtl-Meyer 

angles 

(6) If NPLOT=0, the final output is a printer plot of the blade profile and the flow 
passage. If NPLOT=l, this output is omitted. 
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TABLE n. - SAMPLE OUTPUT 

DESIGN UF SUPERSONIC ILACtS 


CALCULATIONS FGR SUPERSONIC STARTING 
WEI&HT-FLUW PARAMtT ER = C.C5U 

K* ( MA X I = 0. 3447 C= 0.0109 0- 0-4701 H*II(MAXH = 1.8449 

TFE MAXIMUM ObSIGN VALUc FOR VUN) IS 40.5161 UtG WHEN V1LUW) IS 18.0CC0 CEG, V ( CP J IS 59.0000 OEG, GAMMA IS 1.4000 


BETA! IN) = ( 5 . GOO C UtG VlINl = 39.0000 UtG V(UPI 


UfcLTA V = 0.10C0 UEo 


DESIGN PARAMETERS 

= 59.0000 DfcG V (OUT ) * 35.CCCO UEG BETA(OUT) * -68.7395 DEG 

V(ICM) * 18.0000 OEG GAMMA * 1.400C 


LUWEK SURFACE 
UNRCTAIEL TRANSITION ARCS 


INLET N 

X* ( LUW > 

Y* ( LON) 

211 

0. 

0.6725 

201 

— C « C 2 3 / 

L.6727 

151 

— u , G 4 tJC 

0.6 721 

181 

-O.C 726 

C.67I1 

1 ? L 

-0 . C9 6 1 

0. 6695 

161 

-0. 1241 

C. 6675 

151 

-C . 1507 

G . 6650 

141 

-C . 178C 

C..661 9 

131 

-0.2061 

0 • 6 5t)2 

121 

-U .2345 

C. 6539 

111 

-1.2646 

G. 645G 

1 C 1 

— L .2552 

C . 6433 

51 

-0.3268 

L.6 369 

61 

-1.3554 

0.6257 

71 

-0.3531 

G. 621 7 

ei 

- 0 . 42 fcC 

C.6127 

51 

- J . 46 4 3 

C. 602 7 

41 

-C .5015 

G. 5515 

31 

- 0.54 1C 

C. 5 792 

21 

-0.581/ 

C. 5657 

11 

-0.624^ 

L. 5507 

1 

- 0 . 6o86 

C . 5.S41 


Y* (LOW) 

XMLUMI 

OUTLET 

0. 6729 

-0. 

171 

C. 672 7 

0.0237 

161 

0. 6721 

C. 0480 

151 

0.671 L 

0.0728 

141 

0. 6695 

0.0981 

131 

0. 6675 

0. 1241 

121 

0. 6650 

0. 1507 

111 

0.6619 

0.1780 

101 

0.6582 

C.2061 

91 

0.6539 

0.2349 

81 

0.6490 

0.2646 

71 

0. 6433 

0.2952 

fal 

0.6369 

0.3268 

51 

0.6297 

0.3594 

41 

0. 6217 

0.3931 

31 

0.6127 

U . 42 80 

21 

0.6027 

0. 4643 

11 

C. 5915 

0.5019 

1 


KUTATEL ANU TKANSLATtO TRANSITION ARCS 


X* l LOW ) 

- C . 46 < 5 

- C . 4 a A 4 

- C . 5C 1 A 
-C. 5U5 
-C. i J57 
-G . 55 2 G 
-G. 5 7C4 

5c ?s 
-C. tctt 

-C. cj2 3 3 
-0.6412 
-G.t59* 
-L.bl 75 
-C.faGtC 

- C . 7 1 4 1 
-C. M35 
- C . i 52o 
-C. 772C 

- G . / 5 i 5 

- G # fc i l A 
- L « e3 i 5 

- C . « 5 2 G 


Y + ( Loh) 
0. 4841 
0. 4 fa 74 
0. 4302 
0.4322 
G. 4135 
0 . 3 94 C 
0.5 736 
0. 3524 
0. 3303 
u. 3072 
0.2830 
0.257 7 
0.2312 
L.2033 
0. i 741 

t. 1434 

C. 1 11C 
0 . 076 * 
0. 0409 
L .0028 
-0.0375 
- 0 . G 802 


Y*(LOW)-G* 
-0.0703 
-0.0869 
-0.1042 
- 0.1222 
-0.1409 
-0. Ifa04 
-0.1808 
-0.2G19 
— C • <- 24 l 
-j.<l 4 72 
-0.2714 
-0.^967 
-0.3232 
-0.3510 
-O. 58o3 
-0.4110 
-C.4434 
-0.4775 
-0.5135 
-0.5516 
-0.5919 
-0.b34b 


Y»<LUN )- G* 
-C. 1377 
-0.1564 
-C.1758 
-C. I960 
-C.2168 
-C .2385 
-C. 261C 
-L .2843 
-C.iCfct 
-C .33 35 
-C.36C3 
-C.3678 
-C . 4 166 
-C .4466 
-C.47E1 
-G.5U1 
-C. 5457 
-C .5622 


Y* ( LOW ) 
0.4167 
C. 3980 
C.3785 
C . 3584 
C . 3375 
G. 3159 
C.2934 
C.2701 
G . 2458 
L. 2205 
C.1941 
0.1666 
C. 1378 
0.1078 
G.C763 
0.0433 
L • GO 86 
- L. 0278 


X* ( LOW) 
0.52 84 
0.5429 
0.5575 
U.5720 
0.5865 
0.6010 
0.6155 
0.6300 
0.6445 
0 . 6590 
0.6735 
0.6880 
0.7025 
0.71 70 
0.7316 
0. 7461 
0.7607 
0.7753 


OUTLET K 
171 
161 
151 
141 
131 
121 
111 
101 
91 
81 
71 
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TABLE II. - Continued. SAMPLE OUTPUT 


CIhCGLAK akcs 


X*0 ICO*) 

/♦CIIU ) 

- j ,*t6h 

0 . 9o9 1 

-0.9389 

0.9921 

-u.45Qj 

U . 3uC l 

-0.99 1 5 

0 .3J79 

~0.9 

U .3 133 

-d.923b 

0 .3230 

-U.9 1** 3 

0 . 5 30 3 

-L . 9 030 

0.3379 

-0. J933 

U . 5999 

-0. 3dti0 

0 .3512 

-u. 3 

0.5579 

-U.36o5 

0 . 5o9 9 

-u. 3 3t>0 

0.570 7 

-U.3960 

0 .5 766 

— 0 . .1 3 o 3 

0 .5026 

-0. 3262 

0 .3686 

-U. 3i t>9 

0 . 5992 

-0.3033 

0 .5996 

-0.2530 

0 .0096 

-o. 289 9 

U .6099 

-0.2 737 

0 . 0 19b 

-0.2b29 

U .6 194 

-0. 2321 

0 .6239 

-C.291 2 

0 .626 2 

-0.2302 

0.6323 

-G.*191 

0.636 3 

-0.2079 

U .6900 

-0.1967 

U . o9 33 

-U. i 635 

0 . 6*46 9 

-0. 1 792 

0.6500 

— U. 1626 

0 .0529 

-u.i-jh 

U .6557 

-0. 1399 

0 .656 2 

— 0 » l* tts 

0.6666 

—0.1169 

0.662 7 

-0. 1033 

0 .6690 

-O.C93/ 

0 * 666 9 

-0. Cd20 

0.6679 

-0. c 703 

0.6692 

-0. 058o 

0 .670 9 

-o. c909 

U .6713 

-G. 0332 

0 .6 720 

-u.0233 

0.6723 

-0. Cl 1 7 

o .o72d 

w. Gudo 

0 . 6 7 2 9 

0.01 1 7 

0 .6726 

o. 0233 

0.6723 

0. CJ32 

0 .o 72o 

G.09t>9 

0 .o 7 1 3 

u.0386 

O . 6 70 9 

o. 6 /U3 

0 .669 2 

0 . ICIl’O 

0 .66/9 

w.0537 

U .666 9 

o. 1033 

U .6696 

G. 1 169 

0 .062 / 

0. 1/69 

0.6 60 6 

0.1399 

0.6362 

O. 1 319 

0.6337 

0. 16/6 

0 .0 3* 9 

0. 1 79/ 

0 .6300 

0.1 d 35 

U *o9o9 

0.1967 

0 . o 9 3 3 

O.2079 

0.6900 

0.2191 

O .6 36 3 

G.2JU2 

O.o323 

0.2912 

o> *6 2d 2 

0.2321 

0 .6239 

O.2o29 

0.6199 

o.2 7>7 

U .O 19t> 

U. 2 699 

0 .6099 

0.2 950 

0 • 609o 

0.3033 

0.5996 

0.3139 

C : .399 2 

0.3*62 

0.5686 

0.3363 

d .5826 

0 . 3966 

0.5 768 

0.3566 

0 .57o7 

0 . 3 O o 3 

0.369 9 

0.3/63 

0.53/9 

C. JBoo 

u.5312 

0. 3933 

d . 399 9 

0.9030 

0 .3379 

0.9 193 

0 . 5 30 3 

o .4233 

U .3 230 

0.9 320 

0.3153 

C. 99 13 

0 • 5u 7 y 

0 . 9 3o3 

0 . 5ou 1 

0.9369 

0.992 1 

0.9675 

0.9691 

u .9 73d 

0.3738 

U . 9 b 9 i 

0.967b 

0.9921 

0.9369 

o. 3U 01 

0 .43o 3 

0.3079 

0.4915 

0.31 33 

U .4326 

d. 3230 

0.4233 

0.32 69 

0 .4167 


X*U LUk ) 

Y*0 (6Gk J 

-0 . 4 c 7 5 

-O.C703 

-C . 4 5 1 S 

- C • C t 2 * 

-C .9 508 

— 0 . o 5 9 3 

-0.9911 

-C.C965 

-U .9 8*fc 

-C . 0385 

-0.9*36 

-C.C319 

-0.919; 

“ C • 0 * 9 1 

-0.9050 

-0.0170 

-0.3951 

-O.CICL 

~0 .3860 

-C .00 32 

-0.37t3 

0.0C33 

-0 . 3tt 5 

0 .CiCC 

-L . 3 366 

C.olt 3 

-0.3966 

C . C 2 2 9 

-0.3361 

C. C*b9 

-0 « 3 2 6 2 

0.089* 

-0.3155 

0. C 358 

-0.3C85 

C .C952 

-0 . * 5 5C 

0.C1C9 

-0 . 2 E9 9 

0.0515 

-C . 2 7 3 < 

0 .0609 

-0 . * 6 2 5 

0.065 1 

-0.212] 

C . 0 655 

— 0 .2912 

C . 0 7 38 

“0 .2302 

0 . 0 /to 

-C .2 15 1 

C . C8 15 

-0 .20 7S 

0 . 0 8 1 6 

-0.15* 7 

0. C851 

-0 . 1 e 5 5 

0. C5*5 

-0. 1 79* 

0 . C 5 1 6 

“0 .lt*6 

C.C566 

-0 . 1 11 9 

0 . 1C 1 3 

-0 . 1 3 3 S 

0.1038 

-0 * 12 1 9 

0 • 1 0 6 2 

-0 . 1 J t 5 

C . 1 C £ j 

-o . 1 o 1 1 

0.1103 

-C . C 5 J i 

C.1UC 

-U.G8*t 

C. 1135 

— 0 . C 7 C 3 

0.1195 

-t .C 18 t 

0. 116L 

— 0 .C9tS 

0 . 1 its 

-0 . C 3 5 2 

0.11/6 

-0 .0 < 3 5 

0.1181 

-0.0117 

C . 1 i 89 

0 . o C C C 

0 . 1 it 5 

0 . L 1 1 7 

0.1169 

0 .0* J 6 

0 . 1 It 1 

0 .C 882 

C . 1 1 7t 

C.0965 

0 . 1 its 

o .05 8 6 

0.1160 

t.iil: 

0.1195 

0 .C b * 0 

0.1133 

0 .0 5 3 t 

0.1120 

0.1C53 

0. lit 3 

0.1165 

0.1083 

C . 1289 

0 . 1C6* 

C . 133 5 

L. 10 38 

0.1119 

C.K1 j 

0 . ] 62t 

o . 0 5 1 1 

0. 1 

0 . 0 5 It 

0. 1 til 

0 . C S * 5 

0 . 15 6 i 

0 .0 85 1 

0.20 7 5 

0 . <. t It 

0.213] 

0 . C8 IS 

C .* 3C< 

C. 0 i 80 

0 . * 9 1 * 

C . 0 1 3d 

0.25*1 

0. 0 6 S3 

0.2625 

0 .0 6 5 1 

0.2737 

0. OtCS 

0 .*699 

0.0553 

C.2S1L 

O.C 109 

0.3C11 

0 .0 9 5* 

0.3115 

C . 0 3 S e 

0.3*62 

0 .0 892 

O.iitl 

0.0*89 

0 . 3966 

0 .0229 

0.816 6 

0 .0 163 

C .366 1 

0 . C 100 

0.376 3 

C. 0085 

0 . J«it 

- C . 0 C 3 * 

0.3511 

-0.0 ICO 

0 » 4U 5 0 

-o. 01 70 

0.9 193 

- 0 . t * 5 1 

0.9281 

-C.C819 

0.93*6 

- 0 . 0 3 8 5 

0.9911 

-0.0565 

0.95C3 

“0.0 693 

0.9565 

-C.C622 

0.9675 

-C.07C3 

0.9756 

-0.0786 

0 .9t 9 J 

- 0 . C 8 1 S 

0.9521 

- 0 . C 3 5 9 

C.5CC1 

1091 

0 .50 iS 

-0.11*5 

0.5151 

-0.1*18 

u .523C 

-0. 13CS 

0.52t9 

-0.1377 



TABLE H. - Continued. SAMPLE OUTPUT 


UPPER SURFACE 
U NRCT A f E L TRANSITION ARCS 


INLET J 

X M UP ) 

Y* l UP 

201 

0 . 

C. A 8 3 

191 

-t.cito 

0. Aeil 

lei 

-u. C326 

C.AH27 

171 

-C .CAC l 

L . A b2 1 

lei 

- 1 .Cb J 1 

C. Adl2 

151 

-C .C7 7« 

L . A dC C 

1 A 1 

- L .c9 2 1 

C. A7d7 

131 

'■■■ 0 . 1 C 6 C 

C. A7 71 

121 

-t . 1197 

G • 4 753 

111 

-0. 1 3 3 1 

C. A733 

1C1 

-C . 1A62 

C.A711 

91 

- o . i39 1 

C. A6d7 

ei 

-0.171b 

L. Ab62 

71 

- G . 1 6 A 3 

G . A 63 A 

tl 

— 0 « 1 9 b 7 

C. A tL5 

51 

-o .cta9 

C. A3 73 

A 1 

- O . 2 2 G 9 

G. A SAG 

31 

- 0 . 2 3 2 0 

C. A 305 

a 

-0 . 2 A A t 

G » 4 4 6 B 

1 1 

- C . c St i 

C. 4421 9 

l 

- o • 2 ti 7 c 

C.43BB 


Y* (UP) 

X* (UP) 

OUT LEI 

0. 4833 

-0. 

241 

0. 4631 

0.0166 

231 

0. A 82 7 

0.0326 

221 

0. 4821 

0.0481 

211 

0. 4812 

0. 0631 

201 

C. A 800 

0.0778 

191 

0. 4787 

0.0921 

181 

0. 4771 

0. 106 r ' 

171 

0. A 753 

0.1197 

161 

0. A 733 

0.1331 

151 

C.A71 1 

0.1462 

141 

0.4687 

0.1591 

131 

0. 4662 

0.1718 

121 

0. Ab 3 A 

0.1843 

111 

0.4605 

0. 1967 

101 

C. 4573 

0.2089 

91 

0. 4540 

0.2209 

81 

0.4505 

0.2328 

71 

0. AA68 

C. 244b 

61 

0.4429 

0.2563 

51 

0. 4388 

0.2678 

41 

0. 4345 

0.2793 

31 

C.4J01 

0.2906 

21 

0.42 54 

0. 3019 

11 

0. 42C6 

0.3130 

1 


KLi T AT t L AND TRANSLATED TRANSIT IUN ARCS 


*< 


iNLtl J 
<CL 
191 
181 
171 

Lt l 
151 
lAi 
lit 

ia 

in 

1G i 

91 
ai 
71 
6 i 

41 


21 

11 


X* ( UP 1 
-C.34 1 7 

- C . 3 5 l 3 
-C. 3 >4 4 
-G.374S 
-C. 3849 

- L. J94A 
- L. 40 Jt 

- C . A L ^ 3 

-C. 42L / 

- C . 4 2 E t 
-0. AitS 

- C. A A A C 
-l. 4511 

- C. 4 6 tt 

- C . Ati A 7 

- w A 7 i 1 
-C.4772 

- G . 1 d 3 i 

- C. 4389 

- G • 94 4 

- C. A 'AS 7 


Y* luP) 
0* 3 a 1 7 
0. 32 99 
U. 3 l Bi 
U. 3069 
0.2966 
0* 7 BAA 
0.2 7 AA 
u. 2624 
0. 2 3i 6 
0. 24Gb 
U.229B 
U.2 189 
G. 2061 
ij . i 9 n 
0. 1 db5 
u. 1 757 
0. i bA a 
0.1339 
U. 142 9 
0.1319 
C. 12C9 


Y * CUP) ♦G* 
0.89t> 1 
0.6BA3 
v,.«727 
O.bb 13 
0 . o 9 00 
0. dido 
J.82 7d 
0 . B 1 68 
0. 6039 
0.7950 
0. 7d A 1 
*.7733 
C. 7o2 3 
0.7317 
o. 7 A 09 
C. 7 300 
0.7192 
C. 70B3 
0.6973 
0.6B63 
0.6733 


CUP )«G* 

Y* 1 UP) 

C. 697 7 

J. 3A33 

C.EE59 

C.3315 

L • t 7 44 

C. 3200 

C • 8 1 30 

C . 3006 

G . t 5 1 7 

C. 2974 

C .EAC6 

C.2062 

C .8*96 

C.2752 

G . t 1 d 7 

C. 2043 

C . tC 78 

G. 2334 

0.7969 

C. 2426 

0.7861 

C.2317 

C.7753 

C.2210 

C.7646 

C.2102 

C.7538 

G. 1994 

C.7A3C 

C. 1886 

C. 7 322 

0.1778 

C.7213 

0.167C 

C.71C5 

C. 1561 

C • < 5 9 5 

G. 1462 

C.6tt6 

C. 1342 

0.6 7 7 5 

C. 1232 

C ■ 6 6 6 5 

C. 1121 

C.6553 

C. 10G9 

6 • 6 A A 1 

G.C097 

C .6328 

0. G 784 


X*IUP) 

OUTLET 

0.3402 

241 

0.3518 

231 

0. 362 9 

221 

0.3735 

211 

0.3835 

201 

0.3931 

191 

0.4U23 

181 

0.4111 

171 

0.4196 

161 

0.42 77 

151 

0.4355 

141 

0.4430 

131 

0.4502 

121 

0.4571 

ill 

0.4638 

101 

0 . 4 7 03 

91 

0. 4 7b5 

81 

0.4825 

71 

0.4882 

61 

0.4938 

51 

0.4991 

41 

0.5042 

31 

0.5091 

21 

0.5139 

11 

U.51 04 

1 
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TABLE n. - Continued. SAMPLE OUTPUT 

C1HCULAK AKtb 


X*U ( UP) 

Y *0 (UP ) 

X*l< UP J 

Y*CUPHG* 

-0. 34 1 7 

0.341 7 

-0.3417 

0. 646 1 

-U. 3 357 

0 . 3476 

-0.3 38/ 

C.9C2C 

-u.32 l ?b 

0.3634 

-0.3*91 

C . 9 L 7 6 

-0. 32 34 

0 . 3 69 1 

-0.3*24 

C . 4 1 3 6 

-0. 3l 7 l 

0 . 3647 

-0.31? J 

C . 4 l 4 1 

- 0. 3 1 ub 

J . 3 7C 2 

-0. JlUt 

0.4*46 

-0. 3U4 1 

0 . 3766 

-0.3C41 

0.S3C0 

-U.2'775 

u . 3b0 b 

-0.2 4 11 

C .9 25* 

- C. 2 90a 

0 . 3660 

— C . * 40 t 

0 .9403 

-U.2H41 

0.3910 

-0.2841 

0.9454 

-U.2 1 7* 

0.3969 

-0.2772 

0.45C3 

-C.2 tU2 

0 . 4Uu 7 

— 0 .2 ?C* 

C.S 660 

-0.2632 

0.4U63 

-0 .26 3 2 

C. 4697 

-0.2561 

0.4098 

-C .* 6 c 1 

0.9642 

-U.24aV 

0 .4142 

-C.248S 

C . 9 6 1 1 

-O.*4lo 

j • 4 la 6 

-0 .24 It 

C . 9 7 2 9 

-U.234 J 

U .422 7 

-0.2343 
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0 .0834 

1 .0 30 J 
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TABLE H. - Concluded. SAMPLE OUTPUT 



M ISCtLL ANtOJS PARAHtrtHS 

r IHL MINIMUM LUWl K SURFACE PR AKO I L-Mt Yt R ANGLE PKELIcTfcO BY StPARMIUN CUNU111UNS IS IS. 3009 UfcU 

f Hfc MAXIMUM UPPtR SURFACE PR AND J L-Mfc YEP ANGLE PRtOlCTtL BY SEPARATION CONDITIONS IS 56.B766 OfcG 

5 ^ h* ( I N ) * l.e240 MUM = 2.4947 MiUUTI = 2-32*7 MMULD = 1.7b69 

R*(LUrt) = (.6729 M*(LOh» = 1.4tJoO MILCW) = 1.7065 MtUPJ = 3.5297 M* I UP ) = 2.0692 R*(UPl 

j H I TA = 133. 7395 l>EO O* = U.5544 C* * l.fc<c£] 51GMA - 2. 9367 


0 .4B33 
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PROGRAM DESCRIPTION 


Main Program 

The main program generates a table of the inlet and outlet transition arcs of the up- 
per, lower, and translated curves. It also computes several parameters which are per- 
tinent to the blade description, such as Mach numbers and radii, and, as an option, 
plots a blade profile and flow channel. The plotting is done by subroutine PLOTMY, 

(ref. 10). The program variables are 

ALPH fixed angle 

ALPHLN rotation angle for lower- curve inlet transition arc, a f ■ 

ALPHLO rotation angle for lower- curve outlet transition arc a 

l , o 

ALPHUI rotation angle for upper- curve inlet transition arc, a 

U, 1 

ALPHUO rotation angle for upper- curve outlet transition arc, a 

’ u, o 

ALPHUP temporary storage 

ALPLOW temporary storage 

ANGLE logical switch 

BETAN see INPUT 

BETAT outlet flow angle, 

CONVER conversion factor for degrees to radians 

COSALN cosine of ALPHLN 

COSALO cosine of ALPHLO 

COSAUI cosine of ALPHUI 

COSAUO cosine of ALPHUO 

CSTAR blade chord, C* 

DALPH angle increment 

DELF see subroutine ROOT 

DELV see INPUT 

DELXI 1/15 length of straight-line portion of upper-curve inlet arc 
DELXO 1/15 length of straight-line portion of upper-curve outlet arc 
EMJ slope of Mach lines, rm 
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EMK 

EMWJ 

EMWK 

F(V, FN) 

FLO 

FN 

FOFX 

FUP 

GAM 

GAMEXP 

GAMM1 

GAM PI 

GRTY 

GSTAR 

I 

IPRINT 

ISTART 

J 

JJ 

JDEX 

JMAXN 

JMAXO 

JMN 

JN 

JNDEX 

JNN 

JO 

JOO 

K 


slope of Mach lines (eq. (15)), m^ 

slope of wall segments, m. 

slope of wall segments (eq. (17)), m k 

internally defined function (eq. (13)), f(R*) 

see subroutine START 

floating point index 

see subroutine ROOT 

see subroutine START 

ratio of specific heats, y 

1/ (GAM - 1) 

(GAM - l)/2 
(GAM + l)/2 
dummy name 
blade spacing, G* 
counter 
see INPUT 
see INPUT 

variable index for upper curve 
variable index 
variable index 

number of points on upper- curve inlet transition arc 
number of points on upper- curve outlet transition arc 
maximum of JMAXO and JMAXN 
variable index 

number of upper- curve transition arc points to be printed 
variable index 
variable index 
variable index 

variable index for lower curve 


KK 

KKK 

KDEX 

KMAXN 

KMAXO 

KMN 

KN 

KNDEX 

KNN 

KO 

KOO 

KOUNT 

L 

LLL 

LSTORE 

LSTR 

M 

MAXN 

MAXO 

N 

NPER 

NPLOT 

NP1 

NP2 

NP3 

NSUM 

NUM 

P 

PERM 


variable index 

array required by PLOTMY (see ref. 10) 
variable index 

number of points on lower- curve inlet transition arc 
number of points on lower- curve outlet transition arc 
maximum of KMAXN and KMAXO 
variable index 

number of lower- curve transition arc points to be printed 

variable index 

variable index 

variable index 

counter 

counter 

NP1 + NP2 

number of points saved per 5° of turning 

temporary storage 

counter 

integer constant 
integer constant 
variable index 

variable governing selective storage 
see INPUT 

number of points on lower curve which have been saved for plotter 

number of points on translated lower curve which have been saved for plotter 

number of points on upper curve which have been saved for plotter 

total number of points which have been stored for plotter 

counter 

array required by PLOTMY (see ref. 10) 

[(GAM + 1)/(GAM - l)] 1 / 2 
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PHIJ 


PHIJP1 

PHIK 

PHIKP1 

R 

RA 

RECONV 

RIN 

RLOW 

ROUT 

RUP 

SAME 

SIGMA 

SINALN 

SINALO 

SINAUI 

SINAUO 

SM 

SMIN 

SMLOW 

SMOUT 

SMS 

SMUP 

SSMIN 

SSMLOW 

SSMOUT 

SSMUP 

TANBI 

TANBO 


flow direction, < p ^ 
previous value of PHIJ, 
flow direction (eq. (12)), cp ^ 
previous value of PHIK, 

array for storing radii of major vortex- compression- characteristic points, 

array for storing radii of major vortex- expansion- characteristic points, R 

conversion factor for radians to degrees 

1/SSMIN 

jj( 

radius of circular arc of lower curve as calculated in JOKOS, R^ 
1/SSMOUT 

radius of upper- curve circular arc, R* 

see subroutine START 

blade solidity, a 

sine of ALPHLN 

sine of ALPHLO 

sine of ALPHUI 

sine of ALPHUO 

temporary storage for Mach numbers, M 
inlet Mach number, M^ 
lower-surface Mach number, M^ 
outlet Mach number, M Q 
temporary storage for velocity ratio, M* 
upper-surface Mach number, M y 
inlet velocity ratio, M* 
lower-surface velocity ratio, M* 
outlet velocity ratio, M* 
upper-surface velocity ratio, M* 
tangent of BETAN 
tangent of BETAT 
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TEMP 

TEMPP 

TEMPPP 

THETA 

TR 

TX 

TXLO 

TXUP 

TY 

TYLO 

TYUP 

UMJ 

UMJP1 

UMK 

UMKP1 

V 

VIMAX 

VIN 

VLOW 

VLSPMN 

VNL 

VOL 

VOUT 

VUI 


temporary storage 

temporary storage 

temporary storage 

total flow turning angle, 6 

temporary storage for radii 

temporary storage for unrotated x*- coordinates 

temporary storage for values of XLOW(I) 

temporary storage for values of XUP(I) 

temporary storage for unrotated y*- coordinates 

temporary storage for values of YLOW(I) 

temporary storage for values of YUP(I) 

Mach angle, ^ 

previous value of UMJ, jm. +1 

Mach angle (eq. (15)), 

previous value of UMK, M k+1 

temporary storage for Prandtl-Meyer angles, v 

see subroutine START 

see INPUT 

see INPUT 

minimum lower-surface Prandtl-Meyer angle from separation criterion, 

{v l ] ■ 
min 

VIN-VLOW 
VOUT- VLOW 
see INPUT 
VUP-VIN 


VUMAX 

VUP 



see INPUT 
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VUSPMX 

VUT 

XO 

XI 

X2 

XCG 

XCLOW 

XCUP 

XDOWN 

XINTL 

XLOW 

XLOWN 

XLOWO 

XMLOW 

XMUP 

XSIN 

XSOUT 

XUP 

XUPN 

XUPO 

Y AC RSI 
YACRS2 


maximum upper-surface Prandtl-Meyer angle from separation criterion, 
max 

VUP-VOUT 
see subroutine ROOT 
see subroutine ROOT 
see subroutine ROOT 

X*- coordinate of a translated- curve circular arc point 
X*-coordinate of a lower-curve circular arc point 
X*-coordinate of an upper-curve circular arc point 
array for storing X*- coordinate of points to be plotted 
see subroutine ROOT 

array for storing x*-coordinate of unrotated lower transition arc points, 



array for storing X*-coordinate of lower-curve rotated inlet transition arc 

points, (X?) 

L k,i 

array for storing X*-coordinate of lower-curve rotated outlet transition arc 

points, (X*) 

L k, o 

temporary storage for values of -XLOW(I) 
temporary storage for values of -XUP(I) 

X*- coordinate of a point on inlet straight-line portion of upper curve 
X*- coordinate of a point on outlet straight-line portion of upper curve 
array for storing x*-coordinate of unrotated upper transition arc points, 
(x £). 

array for storing X*-coordinate of upper-curve rotated inlet transition arc 
points, (X*) 

j, i 

array for storing X*-coordinate of upper-curve rotated outlet transition arc 
points, (X*) 

j, o 

array for storing Y*- coordinates of points for plotter 
array for temporary storage of Y*-coordinate of translated lower- curve 
points for plotter 
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YCG 

YCLOW 

YCUP 

YLASTI 

YLASTO 

YLOW 

YLOWN 

YLOWO 

YNG 

YSIN 

YSNG 

YSOUT 

YSTG 

YTG 

YUP 

YUPN 

YUPO 


Y*- coordinate of a translated- curve circular arc point 

Y*- coordinate of a lower- curve circular arc point 

Y*-coordinate of an upper-curve circular arc point 

Y -coordinate of first point on inlet side of upper curve 

Y -coordinate of first point on outlet side of upper curve 

array for storing y*- coordinate of unrotated lower transition arc points 

(yf) 

L k 

array for storing Y*- coordinate of lower- curve rotated inlet transition arc 
points, (Yf) 

1 M 

array for storing Y*- coordinate of lower- curve rotated outlet transition arc 
points, (Yf) 

L k, o 

temporary storage for Y*-coordinate of a point on translated- curve inlet 
transition arc 

Y -coordinate of a point on inlet straight-line portion of upper curve 

Y*- coordinate of a point on inlet straight-line portion of translated upper 
curve 

Y -coordinate of a point on outlet straight-line portion of upper curve 

Y -coordinate of a point on outlet straight-line portion of translated upper 
curve 

temporary storage for Y*- coordinate of a point on translated- curve outlet 
transition arc 

array for storing y*- coordinate of unrotated upper transition arc points, (y 

array for storing Y*- coordinate of upper- curve rotated inlet transition arc 
points, (Y*) 

j,i 

array for storing Y*-coordinate of upper-curve rotated outlet transition arc 
points, (Y*) 

h o 


Subroutine ROOT 


Subroutine ROOT is a general routine for finding the roots of equations and is de- 
rived from the "half-interval search" method described in reference 11. This method 
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depends on successively halving an interval which is known to contain the desired root. 
Subroutine ROOT is used to calculate R* from equation (10b). 

A call to ROOT has the form CALL ROOT (X0, X2, XINTL, FOFX, FUNC, XI), 
where the elements of the call vector are 

X0 lower bound of initial root interval 

X2 upper bound of initial root interval 

XINTL initial estimate of value of root 

FOFX given value of dependent variable 

FUNC externally defined function 

XI value of root 

Other program variables are 


A FOFX-F2 

DELF convergence criterion 

F0 function FUNC evaluated at XXO 

F2 function FUNC evaluated at XX2 

FX function FUNC evaluated at X 

KOUNT count of number of iterations performed 

X temporary storage for present estimate of root 

XXO present value of lower bound of root interval 

XX2 present value of upper bound of root interval 


Subroutine START 

Subroutine START is used to compute the maximum value of the inlet Prandtl- Meyer 

angle ( v .) for supersonic starting for various lower- and upper-surface Prandtl- 

1 max 

Meyer angles. Several other parameters of interest are computed by START and are 

printed as output. Among these parameters are the vortex constant for maximum weight 

flow K* , the reduction in weight flow due to two-dimensional flow C, and the maxi- 
max 

mum value of the inlet velocity ratio for supersonic starting (M?) 

max 

A call to START has the form CALL START (VLOW, FLO, VUP, FUP, VIMAX), 
where the elements of the call vector are 

VLOW lower-surface Prandtl-Meyer angle, v ^ 
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FLO 


eq. (10b) evaluated at R* 

VUP upper-surface Prandtl-Meyer angle, 1 / 

FUP eq. (10b) evaluated at R^ 

VIMAX maximum inlet Prandtl-Meyer angle for supersonic starting, (m) 

1 max 

Other program variables are 
BINTGR value of integral (eq. (27b)) 

C reduction in maximum weight flow due to two-dimensional flow (eq. (34b)) 

CINTGR value of integral (eq. (32)) 

FINTL eq. (27) evaluated at K* = XINTL 

max 

F0 eq. (27) evaluated at K* = X0 

max 


FOFX 


eq. (27) evaluated at K* „„ = X2 

max 

eq. (27) evaluated at K* = XAMK 

max 

eq. (34a) 


RATIO ratio of Q to 1 - C 

RLOW radius of circular portion of lower curve, R* 

RUP radius of circular portion of upper curv< , R* 

SAME square of ratio of XAMK to SSMIOW 
SLOPE slope of line 

SSMIAX maximum value of entering velocity ratio for starting, (M*) 

SSMLOW lower-surface velocity ratio, M* 

SSMUP upper- surface velocity ratio, M* 

WSTAR weight-flow parameter (eq. (32)) 

XAMK vortex constant for maximum weight flow, K* 

& ’ max 

XINTL initial estimate of a parameter 
X0 lower bound of a parameter 

X2 upper bound of a parameter 

YINCPT y-intercept of a line 
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Subroutine MSSTAR 


Subroutine MSSTAR is used to determine the minimum lower-surface Prandtl-Meyer 
angle and the maximum upper-surface Prandtl-Meyer angle from separation considera- 
tions. This subroutine uses ROOT and ADSTR in the calculation of these angles. 

A call to MSSTAR lias the form CALL MSSTAR (M, N, VSSTAR), where the ele- 
ments of the call vector are 

M inlet or outlet dimensionless velocity, M* or M*, respectively 

N variable switch 

VSSTAR Prandtl-Meyer angle from separation criterion, (t\) or {v ) 

min max 

The other variables used by MSSTAR are 


A 


B 

C 

D 

MS 

XO 

X2 

XINTL 

FOFX 

SQRDMS 



y - 1 

y + 1 

velocity ratio from separation criterion, (M*) or (M*) 

1 min max 

see subroutine ROOT 
see subroutine ROOT 
see subroutine ROOT 
see subroutine ROOT 
(MS) 2 


Subroutine SIMPS1 


This function subprogram is used to perform numerical integration of explicit func- 
tions of one variable. The integration is performed by a modification of Simpson's rule, 
in which a sufficient number of intervals is used to assure six or more significant fig- 
ures in the result. 

A call to SIMPS 1 has the form ANSWER = SIMPS1 (XMIN, XMAX, FUNC1, KER), 
where the elements of the call vector are 
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XMIN 


lower limit of integration 
XMAX upper limit of integration 

FUNC1 externally defined function of a single variable 
KER storage for flagging result if necessary 
Other program variables are 

A array for storing functional values at certain partition points 

ANS sum of subapproximations 

B array for storing functional values at certain partition points 

C array for storing functional values at certain partition points 

E array for storing difference terms 

FRAC variable tolerance used for subapproximations 

H distance between successive points of partition 

K variable index 

N counter 

NE equivalent to E 

NTEST equivalent to TEST 

P array for storing successive subapproximations 

Q sum of difference terms 

SIMPS 1 value of desired integral 

TEST tester for subapproximations 

T tolerance for difference terms 

V array for storing partition points of interval 


Function Subprograms 

The following function subprograms are used intermittently throughout the main 
program and subroutines: 

FUNCTION ALFUNC (A, B, Y), defined by ALFUNC where 

1 2 ^(y-l) 

ALFUNC = _ (A - BY ; 

Y 
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FUNCTION CFACT (Y), defined by CFACT where 


CFACT = — 
Y 


K 


max 


m: 




FUNCTION QFACT (Y), defined by QFACT where 


QFACT 


« i(r±i . ilzI Y 2 '\ 1/(r " 1) 


Y \ 2 


FUNCTION FRAT (Y), defined by FRAT where 

(y + 1 _ y 


+ 1 y - 1 


FRAT 


2 r 

yy-l ^ y + 1 Y 2 _ y ~ 1 
2 2 


FUNCTION FOFRS (X), defined by FOFRS where 


FOFRS 


^~ + * arc sin ^ — - yj + arc sinj(y + 1)X 2 - yj 


FUNCTION FKMAX (Y, L), defined by FKMAX where 


FKMAX - 



l/(y-l) 


dZ 


afe \ 2 

'M* 


l/(y- 1) 


2 i/(y-i> 

- (1 - Y 2 ) 


and L is an optional switch. 

FUNCTION ADSTR (X), defined by ADSTR where 


ADSTR = 




1 - 




‘-for 


J^X* 


1 + 1 1-£- + 1 


2 \i -2Ulix 2 


y + 1 


(y-l)/y ^ l/2 
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PROGRAM LISTING 


COMMON/ EXP ALF/GAMEXP 
COMMON/ROOTS/DELF 

COMMON/ FAC TOR /PERM , S AWE , GAM , GAMM1 , GAMP 1 , SSMLOW , SSMOP , RECONV , GR TY 

DIMENSION RI800) , R A ( 800 ) , XLOW < 800 > ,YL0W<800) ,XL0WN<800) ,YL0WN(800) 
ltXUP(800) , YUP ( 8 00 ) , XLOWO ( 800 ) , YLOWO ( 800 ) , XUPN ( 800 ) , YUPN ( 800) 

2 XUPO ( 800 ) , YUPO ( 800 ) 

DIMENSION XDOWN ( 400 ) , YACRS1 (400) , YACRS2(200) , KKK ( 1 4 ) , P ( 20 ) 

LOGICAL ANGLE 

EXTERNAL FOFRS 

F(V,FN) = <2.*V) - ( (3. 14159265/2. )* (PERM-1.) ) - ( 2 . * ( FN-1 . > *DEL V ) 

CC INPOT AND TITLE 

I START=0 FOR STARTING AND DESIGN ISTART=1 FOR STARTING ONLY 
N PLOT-O IF PLOT IS DESIRED NPL0T=1 IF PLOT IS NOT DESIRED 
I PR I NT = 0 PRINT ROTATED COORDINATES ONLY IPRINT=1 PRINT ONROTATED 

AND ROTATED COORDINATES 

1 READ (5,11) VIN,VOOT,BETAN,VLOW,VOP,DELV,GAM, I START ,NPLOT , IPRINT 

11 FORMAT ( 7<F6.2,2X),3( II, 2X) ) 

WRITE (6,99) 

99 FORMAT ( 1 H 1 , 38X , 53HD ESIGN OF SUPERSONIC BLA 
IDES) - . 

CC CONVERSION FACTORS AND CONSTANTS 
CONVER = . 174532925E-01 
RECONV = 57.2957796 

C ONE POINT WILL BE PRINTED FOR EVERY NPER POINTS CALCULATED 

IF ( DELV ,GE . 0.2) GO TO 12 
NPER = 10 
GO TO 13 

12 NPER = 1 

13 GAMP 1 = (GAM + 1. )/2. 

GAMM1 = (GAM-1. )/2. 

GAMEXP = 1. /(GAM-1.) 

PERM = SORT ( GAMP1/GAMM1 ) 

DELE = 0.000001 

XO = 1 ./PERM 
X2 = 0.999999999 
XINTL = (XO + X2)/2. 

L S TOR E = ( 5.0/DELV ) /FLOAT ( NPER ) 

DALPH = 1.0*C0NVER 

ANGLE = .TRUE. 

IF ( VLOW .LE. AMIN1 (VUP,VIN,VOUT) ) GO TO 120 
WRITE (6,119) 
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119 FORMAT (//31X,70HV(L0W) MUST BE LESS THAN OR EQUAL TO THE MINIMUM 
1 OF V(UP) , V ( IN) , V ( OUT ) ) 

ANGLE = .FALSE. 


120 IF (VUP .GE. AMAX1 ( VIN,VOUT) ) GO TO 118 
WRITE (6,117) 

FORMAT (//33X,66HV(UP) MUST BE GREATER 


117 


THAN OR EQUAL TO THE MAXIMU 


1M OF V< IN) , V< OUT ) ) 
ANGLE = .FALSE. 


118 VUMAX = ( 3. 16159265/2. )*( PERM-1 . )*RECONV 
IF (VUP , L E . VUMAX) GO TO 116 
WRITE (6,115) VUMAX 

115 FORMAT ( / / A 1 X , 37HV ( UP ) MUST BE LESS THAN V(UPMMAX) 
1 ) 

ANGLE = .FALSE. 


,F9.6,6H DEG 


116 IF (.NOT. ANGLE) GO TO 1 


CC PARAMETERS FOR STARTING 


VLOW = VLOW*CONVER 
FLO = F ( VLOW ,1.0) 

VUP = VUP*CONVER 
FUP = F ( VUP , 1 . 0 ) 

CALL START ( VLOW , FLO , VUP , FUP , V I MAX ) 
IF ( ISTART .NE. 0) GO TO 1 


WRITE (6,97) 

97 FORMAT ( / / 58X , 1 7HDE S I GN PARAMETERS) 

CC MISCELLANEOUS CALCULATIONS 
DEL V = DEL V^CONVt R 
FN = 1 . 

V = V I N*CONVER 
DO A 1=1,2 

FOF X = F(V»FN) 

CALL ROOT ( X0,X2,XINTL,F0FX,F0FRS,X1 ) 

IF (I .EQ. 2) GO TO A 
RIN = XI 

V = VOUT*CONVER 
A CONTINUE 

ROUT = XI 

SSMIN = 1 . /R I N 

CALL MSSTAR ( S SM I N , 0 , VL S PMN ) 

SSMOUT = 1,/ROUT 

CALL MSSTAR ( SSMOUT , 1 ,VUSPMX ) 

SMS = SSMIN 

16 SM~= l SORT<! <l./GAMPn*SMS*SMS)/(l.-CGAMMl/GAMPl)*SMS*SMS) ) 
GO TO ( 17,18,19,20), I 

17 SMIN = SM 
SMS = SSMOUT 


I = 2 
GO TO 16 
18 SMOUT = SM 

TEMP = ( ( ( GAMM1 * SMOUT* SMOUT ) + l. )/( ( GAMM1*SM I N*SM I N ) + 1 . ) )**(GAMP1 
1/(2.*GAMM1) ) 

BETAN = BET AN^CONVER 

BETAT = -ARCOS ( COS ( BETAN ) SM I N/ SMOUT ) *TEMP ) 

BETAT = BETAT #R EC ONV 
BETAN = BETAN*RECONV 
DEL V = DEL V#R ECONV 


CC 


PRINT ALL DESIGN PARAMETERS 
WRITE (6,95) BETAN , V IN , VUP , VOUT , BETAT 
95 FORMAT ( / 2 X , 1 1 HBE T A ( I N ) = ,F7.4,4H DEG , 4X , 8HV ( IN) = 
16X,8HV(UP) = , F8 .4 , AH DEG , 7X , 9H V ( OUT ) = ,F7.4,4H DEG 
2 T ) = ,F8.4,4H DEG) 

WRITE (6,94) DEL V , VLOW, GAM 
94 FORMAT ( / 20X , 1 OHDE L T A V = ,F7.4,4H DEG, 1 1X,9HV ( LOW ) 

1 , 1 1 X , 8HGAMMA = ,F7.4) 


,F7.4,4H DEG, 
,4X, 12HBETA(OU 


= ,F7.4,4H DEG 


CC CONVERT FROM DEGREES TO RADIANS 
VIN = VIN*CONVER 
VOUT = VOUT *C ONVER 
VUP = VUP*CONVER 
VLOW = VLOW*CONVER 
BETAN = BETAN*CONVER 
BETAT = BETAT #CONVER 
DEL V = DEL V*CONVER 


CC CHOOSE LONGEST TRANSITION ARC OF LOWER SURFACE 
VNL = VIN - VLOW 
KMAXN = ( VNL /DEL V ) + 0.5 
VOL = VOUT - VLOW 
KMAXO = ( VOL/DELV ) + 0.5 
KMN = MAXOIKMAXN, KMAXO) 

V = A MAXI ( VIN, VOUT) 

CC CALCULATE R* ( LOW ) =RLOW , M* ( LOW ) =SSMLOW, M(LOW)=SMLOW 
IF (VLOW .EQ. 0.0) GO TO 2 
FN = KMN + 1 
FOF X = F ( V , FN ) 

CALL ROOT (X0,X2,XINTL,F0FX,F0FRS,RL0W> 

GO TO 3 

2 RLOW = 1.0 

3 SSMLOW = l./RLOW 
SMS = SSMLOW 

I = 3 
GO TO 16 
19 SMLOW = SM 


SET INITIAL POINTS FOR LOWER ARC CALCULATIONS 
KNDEX = KMN/NPER 
KDEX = KNDEX 
RA ( KDEX+ 1 ) = RLOW 
XLOW ( KDEX+ 1 ) = 0.0 
YLOW ( KDE X+ 1 ) = RLOW 
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PHIKP1 = -(V-VLOW) + FLOAT (KMN)*DELV 

UMKP1 = ARSIN(SQRT (GAMPl#RLOW#RLOW - GAMM1 ) ) 

TXLO = XL OW ( KDEX+1 ) 

TYLO = YL OW (KDEX+1 ) 

ALPHLN = VNL - BETAN 
ALPHLO = -(VOL+BETAT) 

IF (ALPHLN .LE. 0.0 .AND. ALPHLO ,GE. 0.0) GO TO 110 
ANGLE = .FALSE. 

WRITE (6,111) 

111 FORMAT < //27X, 79HV< LOW ) MUST BE GREATER THAN OR EQUAL TO V ( I N ) - B 
1ETA(IN) AND V(OUT) + BET A ( OUT ) ) 

CC CHOOSF LONGEST TRANSITION ARC OF UPPER SURFACE 
110 VUT = VUP - VOUT 

JMAXO = ( VUT/DELV )+0.5 
VUI = VUP - VIN 
JMAXN = ( VU I / DEL V ) +0 . 5 

JMN = MAXO ( JMAXO, JMAXN ) 

V = AMIN1 (VOUT, VIN) 

CC CALCULATE R* ( UP ) =RUP , M* ( UP ) = SSMUP , M(UP)=MUP 
FN = - ( J MN+ 1 ) + 2 
FOF X = F ( V , FN ) 

CALL ROOT <X0,X2,XINTL,F0FX,F0FRS,RUP) 

SSMUP = l./RUP 
SMS = SSMUP 
I = A 
GO TO 16 
20 SMUP = SM 

CC SET INITIAL POINTS FOR UPPER ARC CALCULATIONS 
JNDEX = JMN/NPFR 
JDEX = JNDEX 
R( JDEX+1 ) = RUP 

XUP< JDEX+1 ) = 0.0 
YUPIJDEX+1) = RUP 

PHIJP1 = - ( VUP- V ) + FLOAT ( JMN ) *DELV 
UMJP1 = ARS I N ( SORT ( GAMP1*RUP*RUP - GAMM1 ) ) 

TXUP = XUP (JDEX+1 ) 

TYUP = YUP (JDEX+1 ) 

ALPHUI = VUI - BETAN 
ALP HUO = - ( VUT + BETAT ) 

IF (ALPHUI .LE. 0.0 .AND. ALPHUO .GE. 0.0) GO TO 112 
ANGLE = .FALSE. 

WRITE (6,113) 

113 FORMAT ( / / 2 8X , 7 5HV ( UP ) MUST BE LESS THAN OR EQUAL TO V < I N ) + BET A ( 
1 1 N ) AND V ( OUT ) - BE T A ( OUT ) ) 

112 IF (.NOT. ANGLE) GO TO 1 

IF (VIN .EG. VL OW .AND. VLOW .EQ. VOUT) GO TO 100 

(^♦♦♦CALCULATE COORDINATES FOR LOWER TRANSITION ARC - UNROTATED 
K.OEX = K.NDEX + 1 
NUM = 0 
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V = AMAX1 ( VINtVOUT ) 

DO 30 KK = 1 * KMN 
K = (KMN+1 ) - KK 
NUM = NUM + 1 
PHIK = PHIKP1 - DEL V 
FN = K 

FOF X = F(VtFN) 

CALL ROOT (X0,X2,XINTL,F0FX,F0FRS,TR) 

TX = TR*SIN(PHIK) 

TY = TR*COS(PHIK) 

EMWK = T A N ( —PH I K P 1 ) 

UMK = ARSIN(SORT (GAMP1*TR*TR - GAMM1 ) ) 

EMK = -TAN( ( PHIK+UMK+PHIKP1+UMKP1 ) /2. ) 

TEMP = TYLO - EMWK*TXLO 

TEMPP = TY - 6MK*TX 

TEMPPP = EMK - EMWK 

TXLO = (TEMP - TEMPP ) /TEMPPP 

TYLO = <(EMK*TEMP) - ( EMWK*TEMPP ))/ TEMPPP 

PHIKP1 = PHIK 

UMK P 1 = UMK 

CC SAVE EVERY =NPER-T H= POINT 

N = NUM - (NUM/NPER ) *NPER 
IF (N .GT. 0) GO TO 30 
KDEX = KDEX - 1 
RA(KDEX) = TR 
XLOW(KDEX) = TXLO 
YLOW(KDEX) = TYLO 

30 CONTINUE 

C****CALCULATE COORDINATES FOR LOWER TRANSITION ARC - ROTATED 
100 KDEX = KNDEX + I 
KMN = KMN + I 
SINALN = SIN(ALPHLN) 

COS AL N = COS(ALPHLN) 

SINALO = SIN(ALPHLO) 

COSALO * COS(ALPHLO) 

KN = (KMAXN/NPER ) + 2 
KO = ( KMA XO/NPER ) + 2 

DO 40 KK=I,KDEX 
K = (KDEX+1) - KK 
KN = KN - I 
KO = KO - 1 

IF (KN .LE. 0) GO TO 401 

XLOWN(KN) = YLOW(K )*SINALN + XLOW ( K ) *COSALN 
YL OWN ( K N ) = YLOW(K)*COSALN - XLOW ( K ) *S I NALN 
401 IF (KO .LE. 0) GO TO 40 

XLOWO(KO) = YLOW(K)*SINALO - X LOW ( K ) *COS ALO 
YLOWO(KO) = YLOW(K)*COSALO + XL OW ( K ) *S I NALO 
40 CONTINUE 

IF (VIN .EO. VUP .AND. VUP .EO. VOUT ) GO TO 200 
C****CALCULATE COORDINATES FOR UPPER TRANSIIION ARC - UNROT ATE D 
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JDEX = JNDEX + 1 
NUM = 0 

V = A M I N 1 ( VOUT ♦ V I N ) 

DO 41 JJ=1 ♦ JMN 

J = (JMN+1) - JJ 

NUM = NUM + 1 

PHIJ = PHIJP1 - DEL V 

FN = -J + 2 

EOF X = F(VtPN) 

CALL ROOT (XO t X2,XINTL,FOFX f FOFRS,TR) 

TX = TR*SIN(PHIJ) 

TY = TR*C0S(PH1J) 

EMWJ = TAN(-PHIJPl) 

UMJ = ARSIN(SQRT(GAMP1*TR*TR - GAMM1 ) ) 

EMJ = TAN( <-PHIJ+UMJ-PHIJPl+UMJPl>/2.) 

TEMP = TYUP - EMWJ*TXUP 

TEMPP = TY - EMJ*TX 

TEMPPP = EMJ - EMWJ 

TXUP = (TEMP - TEMPP) /TEMPPP 

TYUP = MEMJ*TEMP) - ( EMWJ*TEMPP ))/ TEMPPP 

PHIJP1 = PHIJ 

UMJP1 = UMJ 

CC SAVE EVERY =NPER-TH= POINT 

N = NUM - ( NUM/NPER ) *NPE R 
IF (N .GT. 0) GO TO 41 
JDEX = JDEX - 1 
R(JDEX) = TR 
XU P (JDEX ) = TXUP 
YUP (JDEX) = TYUP 
41 CONTINUE 

C****CALCULATE COORDINATES FOR UPPER TRANSITION ARC - ROTATED 
200 JDEX = JNDEX + 1 
JMN = JMN + 1 
SINAUI = SIN( ALPHUI ) 

COS AU I = COS ( ALPHUI ) 

SINAUO = SIN(ALPHUO) 

COSAUO = COS(ALPHUO) 

JN = ( J MAXN/NPER ) + 2 
JO = (JMAXO/NPER) + 2 

DO 47 JJ=1 , JDEX 
J = ( JDEX+1 > - JJ 

JO = JO - 1 
JN = JN - 1 

IF (JO .LE. 0) GO TO 471 

XUPO(JO) = YUP ( J >*SINAUO - XUP ( J ) *COSAUO 
YUPO(JO) = YUP ( J ) *COSAUO + XUP ( J ) *S I NAIJO 
471 IF (JN * L E . 0) GO TO 47 

XUPN(JN) = YUP< J )*SINAIJI + XUP ( J ) *C0 S AU I 
YUPN(JN) = YUP ( J ) *COSAU I - XUP ( J ) *S I NAU I 
47 CONTINUE 

CC CALCULATE G* - THE DIMENSIONLESS BLADE SPACING 
TANBI = TAN(BETAN) 

YL AST I = YUPN(l) + T ANB I * ( X LOWN ( 1 ) - XUPN(l)) 
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GST AR 


YLOWN ( 1 ) 


YLAST I 


CC TITLES 

WRITE (6,93) 

93 FORMAT < / / / 56X , 2 5FIL 0 W E R SURFACE) 

IF ( IPRINT .EQ. 0) GO TO 866 
WRITE (6,88) 

88 FORMAT ( / 66 X , 2 5HUNR0 T A T E D TRANSITION ARCS) 

WRITE (6,87) 

87 FORMAT (//2X,8H INLET K ,6X , 8HX* ( LOW ) , 3X , 8HY* ( LOW ) , 68X , SHY* ( LOW ) 

I , 3X , 8HX* ( LOW ) , 2 X , 8H0UTLET K) 

C****PRINT COORDINATES FOR LOWER TRANSITION ARC - UNROTATED 
KDEX = KNDEX + 2 
DO 51 KK=1 ,KMN,NPER 
K = (KMN+1 ) - KK 
KN = (KMAXN+2) - KK 
KO = ( KMAXO+2 ) - KK 
KDEX = KDEX - 1 
XMLOW = —XL OW (KDEX) 

IF (KN .GT. 0 .AND. KO .GT. 0) GO TO 510 

IF (KN .LE. 0) GO TO 511 

WRITE (6,861) KNfXLOW(KDEX) » YLOW ( KDEX ) 

861 FORMAT ( 6 X , I 6 , 5 X , F 8 . 6 , 3 X , F 8 . 6 ) 

GO TO 51 

511 WRITE (6,863) YL OW ( K DE X ) , XML OW , KO 
863 FORMAT ( 1 00 X , F 8 . 6 , 3 X , F 8 . 6 , 5 X , 1 6 ) 

GO TO 51 

510 WRITE (6,860) KN , XLOW ( K DE X ), YLOW ( KDEX ), YLOW I KDE X ), XMLOW , KO 
860 FORMAT ( 6 X , I 6 , 5 X , F8 . 6 , 3 X , F 8 . 6 , 6 8 X , F 8 . 6 , 3 X , F8 . 6 , 5 X , I 6 ) 

51 CONTINUE 

CC TITLES 

866 WRITE (6,92) 

92 FORMAT ( 1 HL , 66X , 38HR0T AT ED AND TRANSLATED TRANSITION ARCS) 

WRITE (6,86) 

86 FORMAT (//2X,8H INLET K , 6X , 7HX* ( LOW ) , 7 X , 7H Y* ( LOW ) , 6 X , 10H Y* ( LOW ) -G* 
1,28X, 10HY*(L0W)-G*,5X,7HY*( LOW ) , 7X , 7H X* ( LOW ) , 5X , 8H0UTL E T K) 

C****PRINT COORDINATES FOR LOWER TRANSITION ARC - ROTATED 
M = 1 

XDOWN(M) = XLOWN(l) 

YACRS1 ( M) = YLOWN ( 1 ) 

YACRS2 ( M ) = YLOWN ( 1 ) - GSTAR 

M = M + 1 

X DOWN ( M ) = XLOWO(l) 

YACRSl(M) = YLOWO(l) 

YACRS2 ( M ) = YLOWO(l) - GSTAR 

CC STORE POINTS FOR PLOTTER - ONE POINT FOR EVERY FIVE DEGREES OF TURNING 
MAXN = ( KMAXN/NPER ) + 1 
MAXO = ( KMAXO/NPER ) + 1 
KNN = MAXN + 1 
KOO = MAXO + 1 
I = 0 


48 



DO 55 KK=1 ,KMN,NPER 
KN = (KMAXN+2) - KK 
KO = (KMAXO+2) - KK 
KNN = KNN - 1 
KOO = KOG - 1 
1 = 1 + 1 
LSTR = LSTORE*I 

IF (KN .GT. 0 .AND. KO .GT. 0) GO TO 550 
IF (KN .LE. 0) GO TO 551 
IF (LSTR .GT. MAXN) GO TO 559 
M = M + 1 

X DOWN ( M ) = XLOWN ( LSTR ) 

YACRS1 ( M ) = YLOWN(LSTR) 

YACRS2 ( M ) = Y LOWN (LSTR) - GSTAR 

559 YNG = YLOWN ( KNN ) - GSTAR 

WRITE (6,831) KN, XLOWN ( KNN ), YLOWN ( KNN ), YNG 
831 FORMAT ( AX , I 4, 7X , E8 . A , 6 X , F 8 .4 , 6X , F8 . A ) 

GO TO 55 

551 IF (LSTR .GT. MAXO) GO TO 557 
M = M + 1 

X DOWN ( M ) = XLOWO (LSTR) 

YACRSl(M) = YLOWO (LSTR) 

YACRS2 ( M ) = YLOWO (LSTR) - GSTAR 

557 YTG = YLOWO(KOO) - GSTAR 

WRITE (6,833) YTG , YLOWO ( KOO ), XLOWO ( KOO ) ,K0 
833 FORMAT ( 8 1 X , F8 . A , 6 X , F 8 . A , 6 X , F 8 . A , 7 X , I A ) 

GO TO 55 

550 YNG = YLOWN ( KNN ) - GSTAR 
YTG = YLOWO ( KOO ) - GSTAR 
IF (LSTR .GT. MAXN) GO TO 558 
M = M + 1 

XDOWN(M) = XL OWN (LSTR) 

YACRSl(M) = YLOWN (LSTR) 

YACRS2 ( M ) = YLOWN ( LSTR ) - GSTAR 

558 IF (LSTR .GT. MAXO) GO TO 556 
M = M + 1 

XDOWN(M) = XLOWO (LSTR) 

YACRSl(M) = YLOWO(LSTR) 

Y ACRS2 ( M ) = YLOWO ( LSTR ) - GSTAR 
556 WRITE (6,830) KN , XLOWN ( KNN ), YLOWN ( KNN ), YNG , Y TG, YLOWO ( KOO ) 

830 FORMAT ( AX , I A, 7 X , F8 . A , 6 X , F8 . A , 6X , F8 . A , 30X , F8 . A , 6X , F 8 . A , 6X , 
1 1 A ) 

55 CONTINUE 
M = M + 1 

XDOWN(M) = XLOWN (MAXN ) 

YACRSl(M) = YLOWN ( MAXN ) 

YACRS2 < M ) = YLOWN ( MAXN ) - GSTAR 

M = M + 1 

XDOWN(M) = XLOWO ( M AXO ) 

YACRSl(M) = YLOWO ( MAXO ) 

YACRS2 ( M ) = YLOWO ( MAXO ) - GSTAR 

q****cirCULAR ARC (LOWER) 

IF ( I PR I NT .EQ. 0) GO TO 810 


XLOWO ( KO 
F8.4,7X, 
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WRITE (6,81) 

81 FORMAT ( //60X , 1 3HC I RCUL AR ARCS/ /60X,8HX*C ( LOW ) , 3X, 8H Y*C ( LOW ) , 16X, 
1 8HX*C(LOW) ,3X,11HY#C(L0W)-G*) 

810 M = M + 1 

X DOWN ( M ) = 0.0 

YACRSl(M) = R LOW 

YAC RS2 ( M ) = RLOW - GSTAR 

THETA = (BETAN - BET AT ) *RECONV 

ALPH = ALPHLO + DALPH 

ALPLOW = ALPHLN 

KOUNT = 0 

60 XCLOW = RLOW*SIN( ALPLOW ) 

YCLOW = RLOW*COS ( ALPLOW ) 

XCG = XCLOW 

YCG = YCLOW - GSTAR 

KOUNT = KOUNT + 1 

IF (KOUNT .NE. LSTORE ) GO TO 601 

KOUNT = 0 

M = M + 1 

X DOWN ( M ) = XCLOW 

YACRSl(M) = YCLOW 

YACRS2 ( M ) = YCG 

601 IF (IPRINT .EQ. 0) GO TO 800 

WRITE (6,80) XCLOW, YCLOW, XCLOW, YCG 

80 FORMAT ( 39 X , F 8 . 6 , 3 X , F 8 .6 , 1 6X , F 8 . 6 , 3X , F8 . 6 ) 

800 ALPLOW = ALPLOW + DALPH 

IF ( ABS ( ALPH-ALPLOW) .LF. 0.001) GO TO 56 

IF (ALPHLO .LT. ALPLOW .AND. ALPLOW .LT. ALPH) ALPLOW = ALPHLO 
GO T 0 60 

CC STORE THE TRANSLATED LOWER ARC FOR PLOTTER 

56 NP1 = M 

DO 2000 1=1, NP1 
M = M + 1 

X DOWN ( M ) = XDOWN ( I ) 

2000 YACRSl(M) = YACRS2U) 

NP2 = NP1 

CC TITLES 

WRITE (6,79) 

79 FORMAT ( 1 H 1 , 53X , 2 5HU PPER SURFACE) 

IF (IPRINT .EQ. 0) GO TO 700 
WRITE (6,76) 

76 FORMAT ( / 56X , 25HUNR0T AT ED TRANSITION ARCS) 

WRITE (6,73) 

73 FORMAT (//2X,8H INLET J,3X,8H X*(UP),3X,8H Y* ( UP ) , 68 X , 8H Y* ( UP ) 

1 , 3X , 8H X*(UP) ,3X,8H0UTLET J) 

C****pr I NT COORDINATES FOR UPPER TRANSITION ARC - UNROTATED 
JDEX = JNDEX + 2 
DO 65 JJ=1 ,JMN,NPER 
J = (JMN+1) - JJ 
JO = (JMAXO+2) - JJ 
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JN = (JMAXN+2) - JJ 
JDE X = JDEX - 1 
XMUP = -XUP(JOEX) 

IF (JN .GT. 0 .AND. JO .GT. 01 GO TO 650 

IF (JN .IE. 0) GO TO 651 

WRITE (6,720) J N , XUP ( J DE X ) , YUP ( JDE X ) 

720 FORMAT ( 4X , I4,5X,F8.4,3X»F8.4) 

GO TO 65 

651 WRITE (6,723) YUP ( J DE X ) , XMUP , JO 
723 FORMAT ( 1 OOX , F 8 .4 , 3X , F 8 . A , 5X , I 4 ) 

GO TO 65 

650 WRITE (6,721) JN , XUP ( J DE X ) , YUP ( J DE X ) , YUP ( J DE X ) , XMUP , JO 

721 FORMAT ( AX , I 4 , 5X , F8 . 4 , 3 X , F 8 .4 , 68 X , F 8 . 4 , 3X , F 8 .4 , 5X , I 4 ) 

65 CONTINUE 

CC TITLFS 

700 WRITE (6,78) 

78 FORMAT ( 1HL , 46X , 38HROT ATED AND TRANSLATED TRANSITION ARCS) 

WRITE (6,70) , _ „„ 

70 FORMAT <//2X,8H INLET J,6X,7H X*(UP),7X,7H Y*(UP),6X, 9HY* ( UP ) +G* , 
129X,9HY*(UP)+G*,6X,7H Y*(UP),7X,7H X*(UP)»5X, 8H01JTL ET J) 

L = NP1 + NP2 

c ****PR1NT COORDINATES FOR UPPER TRANSITION ARC - ROTATED 

CC STORE POINTS FOR PLOTTER - ONE POINT FOR EVERY FIVE DEGREES OF TURNING 
L = L + 1 

XDOWN(L) = XUPN(l) 

YACRSl(L) = YlJPN(l) 

L = L + 1 

XDOWN(L) = XUPO< 1 ) 

YACRSl(L) = YUPO(l) 

MAXO = ( JMAXO/NPER ) + 1 

MAXN = ( JMAXN/NPER ) + 1 
JOO = MAXO + 1 
JNN = MAXN + 1 
I = 0 

DO 303 JJ = 1 , JMNfNPER 
JO = ( JMAXO+2 ) - JJ 
JN = (JMAXN+2) - JJ 
JOO = JOO - 1 
JNN = JNN - 1 
1 = 1 + 1 
LSTR = L S TORE* I 

IF (JN .GT. 0 .AND. JO .GT. 0) GO TO 3030 
IF (JN .LE. 0) GO TO 3031 
IF (LSTR .GT. MAXN) GO TO 688 
L = L + 1 

XDOWN(L) = XUPN(LSTR) 

YACRSl(L) = YUPN(LSTR) 

688 YNG = YUPN(JNN) + GST AR 

WRITE (6,68) JN,XUPN( JNN) ,YUPN( JNN) ,YNG 
68 FORMAT ( 4X , I 4 , 7 X , F8 . 4 , 6 X , F 8 . 4 , 6 X , F 8 . 4 ) 

GO TO 303 

3031 IF (LSTR .GT. MAXO) GO TO 689 


51 



L = L + 1 

XOOWN(L) = XUPO(LSTR) 

YACRSl(L) = YUPO(LSTR) 

689 YTG = YUPO(JOO) + GSTAR 

WRITE (6,683) Y TG , YIJPO ( JOO ) , XUPO ( JOO ) , JO 
683 FORMAT ( 8 1 X , F 8 . 4 , 6 X , F 8 . 4 , 6 X , F 8 . 4 , 7X , I 4 > 

GO TO 303 

3030 YNG = YUPN(JNN) + GSTAR 
YTG = YOPO(JOO) + GSTAR 
IF (LSTR .GT. MAXN) GO TO 670 
L = L + 1 

XDOWN(L) = XUPN(LSTR) 

YACRSl(L) = YUPN(LSTR) 

670 IF (LSTR .GT. MAXO) GO TO 671 
L = L + 1 

XDOWN(L) = XUPO (LSTR) 

YACRS1 ( L ) = YUPO( LSTR ) 

671 WRITE (6,680) J N , XUPN ( JNN ) , YUPN ( J NN ) , YNG , YTG , 
1 JO 

680 FORMAT ( 4X , I 4 , 7 X , F8 . 4 , 6 X , F 8 . 4 , 6 X , F8 . 4 , 30X , F8 . 
114) 

303 CONTINUE 
L = L + 1 

XDOWN(L) = XUPN(MAXN) 

YACRSl(L) = YUPN(MAXN) 

L = L + 1 

XDOWN(L) = XUPO ( MAXO ) 

YACRSl(L) = YUPO(MAXO) 

C****CIRCULAR ARC (UPPER) 

IF (IPRINT .EG. 0) GO TO 6700 
WRITE (6,67) 

67 FORMAT ( / / 60X, 13HCIRCULAR ARCS//40X , 8HX*C ( UP ) 
1 8HX*C(UP) ,3X,10HY*C(UP)+G*) 

6700 L = L + 1 

XDOWN(L) = 0.0 
YACRSl(L) = RUP 

ALPH = ALPHUO + DALPH 
ALPHUP = ALPHUI 
KOUNT = 0 

305 XCUP = RUP*SIN< ALPHUP) 

YCUP = RUP*COS( ALPHUP) 

XCG = XCUP 

YCG = YCUP + GSTAR 

KOUNT = KOUNT + 1 

IF (KOUNT .NE. L STORE ) GO TO 672 
KOUNT = 0 
L = L + 1 
XDOWN(L) = XCUP 
YACRSl(L) = YCUP 

672 IF (IPRINT .EG. 0) GO TO 660 


YUPO ( JOO) , XUPO (JOO) , 
4,6X,F8.4,6X,F8.4,7X, 


,3X,8HY*C(UP) , 1 6 X , 


52 



WRITE (6,66) XCUP, YCIJP, XCUP, YCG 
66 FORMAT ( 39 X , F 8 . A , 3 X , F 8 . 6 , 1 6 X , F 8 . A , 3X , F8 . 6 ) 

660 ALPHUP = ALPHUP + DALPH 

IF (ABS< ALPH-ALPHUP) .LF. 0.001) GO TO 306 

IF (ALPHUO . L T . ALPHUP .AND. ALPHUP . L T . ALPH) ALPHUP = ALPHUO 
GO TO 305 

C****CALCULATE COORDINATES FOR STRAIGHT LINE PORTION OF UPPER ARC 
CC FIFTEEN POINTS ARE CALCULATED FOR PLOTTING PURPOSES 

306 IF ( I PR I NT .EG. 0) GO TO 3070 
WRITE (6,307) 

307 FORMAT ( / / 59X , 1 AHS TR A I GHT LINES//5X,8H X*S ( IN) ,5X,8H Y*S(IN) ,3X, 
1 10HY*S ( I N ) +G* , 5AX , 1 1HY*S (OUT ) +G* , 2 X , 8HY*S ( OUT ) , 5X , 8HX*S ( OUT ) ) 

3070 KOUNT = -1 

DELXI = ( XUPN(l) - XLOWN(l) )/15. 

DEL XO = ( XLOWO(l) - XUPO(l) )/)5. 

XSIN = XUPN ( 1 ) 

YSIN = YUPN ( 1 ) 

XSOUT = XUPOI 1 ) 

YSOUT = YUPO( 1 ) 

TANBO = TAN(BETAT) 

GO TO 309 

310 XSIN = XSIN - DELXI 
XSOUT = XSOUT + DELXO 

YSIN = YUPN ( 1 ) + T ANB I *( XSIN - XUPN ( 1 ) ) 

YSOUT = YUPO(l) + TANBO*< XSOUT - XUPO(l) ) 

309 YSNG = YSIN + GST AR 
YSTG = YSOUT + GST AR 
IF (XSIN .LE. XLOWN(l) ) GO TO 312 
KOUNT = KOUNT + 1 
N = KOUNT - ( KOUNT / 3 ) *3 
IF (N .GT. 0) GO TO 673 
L = L + 1 
XDOWN(L) = XSIN 
YACRS1 (L) = YSIN 

673 IF (IPRINT .EG. 0) GO TO 3133 
WRITE (6,313) XSIN, YSIN, YSNG 
313 FORMAT ( 5X , F8 . A , AX , F 8 . A , AX , F 8 . A ) 

3133 IF (XSOUT .GE. XLOWO(l) ) GO TO 310 
IF (N .GT. 0) GO TO 67A 
L = L + 1 
XDOWN(L) = XSOUT 
YACRSl(L) = YSOUT 
67A IF (IPRINT .EG. 0) GO TO 310 

WRITE (6,315) YSTG, YSOUT, XSOUT 
315 FORMAT ( 1 H+ , 9 3 X , F 8 . A , A X , F 8 . A , A X , F 8 . A ) 

GO TO 310 

312 IF (XSOUT .GE. XLOWO(l) ) GO TO 311 
IF ( IPRINT .EG. 0) GO TO 310 
WRITE (6,321) YSTG, YSOUT, XSOUT 
321 FORMAT ( 9 AX , F8 . A , A X , F 8 . A , AX , F 8 . A ) 

GO TO 310 

311 NP 3 = L - (NP1 + NP2 ) 


NSUM = NP1 + NP2 + NP3 + 1 


XDOWN(NSUM) = 0.0 
YACRSl(NSUM) = 0.0 

C****MI SCELLANEOUS C ALCULAT IONS 
WRITE (6,622) 

622 FORMAT (//54X,24HMISCELLANE0US PARAMETERS//) 

YLASTO = YUPO(l) + TANBO* ( XLOWO ( 1 ) - XUPO(l)) 

CSTAR = SQR T ( ( < XLOWO ( I ) - XLOWN ( 1 ) ) **2 > + ((YLOWO(l) - YLOWN(l)) 

1**2 ) ) 

SIGMA = CSTAR /GSTAR 
WRITE (6,999) VLSPMN, VUSPMX 

999 FORMAT <17X,84HTHE MINIMUM LOWER SURFACE PRANOTL-MEYER ANGLF PREDI 
1CTED BY SEPARATION CONDITIONS IS ,F9.4,4H DEG/ / 1 7X , 84HTHE MAXIMUM 
2UPPER SURFACE PRANDTL-MEYER ANGLE PREDICTED BY SEPARATION CONDITIO 
3NS IS , F 9 . 4 , 4H DEG) 

WRITE (6,1000) SSMIN, SM I N , SMOUT , S SMOUT 

1000 FORMAT (/25X, 9HM*(IN) = ,F8.4,2X,9H M ( I N ) = , F8 .4 , 10X , 9HM ( OUT ) = 
1,F8.4,5X,10HM*(0UT ) = ,F8.4) 

WRITE (6,1001) R LOW ,SSMLOW,SMLOW,SMUP, SSMUP , RUP 

1001 FORMAT ( / 2 X , 9EIR* ( LOW ) = , F 8 . 4 , 5 X , 1 OHM* ( LOW ) = , F 8 . 4 , 2 X , 9HM ( LOW ) = , 

1F8.4,10X,9H M ( UP ) = ,F8.4,5X,10H M*(UP) = , F 8 . 4 , 2 X , 9HR* ( UP ) = , 
2F8.4 ) 

WRITE (6,1002) THETA, GSTAR, CSTAR, SIGMA 

1002 FORMAT ( / 1 1 X, 8HTHETA = ,F8.4,4H DEG,12X,5HG* = , F8 . 4, 1 3X , 5HC* = , 
1F8.4,11X,8HSIGMA = ,F8.4) 

IF (NPLOT .NE. 0) GO TO 1 

CC IF PLOTMY IS NOT AVAILABLE, REMUVE THE FOLLOWING CARDS 
C****MULT I PLE PLOT - START 
LLL = NP1 + NP 2 

CALL SORTXY ( XDOWN ( 1 ) , YACRS 1 ( 1 ) , NP 1 > 

CALL SORTXY ( XDOWN ( NP 1+ 1 ), YAC RS 1 ( NP 1+ 1 ) ,NP2 ) 

CALL SORTXY ( XDOWN ( LL L+ 1 ), YACRS 1 < L LL+ 1 ), NP3 ) 

P ( 1 ) = 5.0 
P ( 3 ) = 12.0 
P ( 4 > = 20.0 

P < 1 1 ) = ( ( 1 . - AMIN1 ( YACRS1 ( 1 ) ,YACRS1 (NP1 + 1 ) ,YACRS1 ( LLL + 1 ) ) ) /100. )* 
1 ( 10 .**4) 

P ( 6 ) = 2.0 

P ( 7 ) = AMIN1( XDOWN( 1 ) ,XD0WN(NP1 + 1 ) , XDOWN (LLL+1 ) ) *(10.**4) 

P ( 8 ) = P( 11 )*( 5./3. ) 

P ( 9 ) = 2.0 

P(10) = AMIN1 ( YACRS1 ( 1 ) ,YACRS1 (NP1+1 ) ,YACRS1 ( LLL+1 ))*( 10. **4) 

KKK ( 1 ) = 55 
KKK< 2 ) = 4 

KKK ( 3 ) = NP 1 
KKK ( 5 ) = NP2 
KKK ( 7 ) = NP 3 

KKK ( 9 ) = 1 

DATA KKK (4) ,KKK(6),KKK<8 ) / 1 H* , 1 H* , 1 H+/ , KKK ( 10 > / 1 HO / 

CALL PLOTMY ( XDOWN , YAC R S 1 , KKK , P ) 

C****MUL T I PL E PLOT - STOP 
GO TO 1 

END 
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$ I BF TC. ROO 


LIST 


C 

c 

c 

c 

c 

c 


SUBROUTINE ROOT ( XO , X2 » X I NT L , FOF X , FUNC , X 1 ) 

COMMON/ ROOTS/DELF 

DOUBLE PRECISION X,XXO,XX2 


X SUCH THAT FUNC(X) = FOFX WHERE FOFX IS A KNOWN 


OR ( F X » F2 ) WHERE FX IS THE PREVIOUS 


WE ARE SEEKING AN 
FUNCTIONAL VALUE 

1 LOCATE FOFX IN ( FO , F X > 

APPROXIMATION TO FOFX 

2 LET X = 1 / 2 ( X XO+ X ) OR X = 1/2(X+XX2) 

3 IS FUNC ( X ) = FOFX = IF NOT, REPEAT PROCEDURE 


XXO = XO 
XX2 = X2 
FO = FUNC (XXO) 

F2 = FUNC ( X X2 ) 

IF ( FOFX .LT. FO .AND. FOFX .LT. F2 .OR. FOFX .GT. FO .AND. 
1 FOF X .GT. F2 ) GO TO 1005 

IF ( ABS ( FOFX-FO ) .LE. DELF) GO TO 1007 
IF ( ABS ( F0FX-F2 ) .LF. DELF) GO TO 1008 


X = XINTL 
KOUNT = 0 

1000 XI = X 

KOUNT = KOUNT + 1 
A = FOFX - F2 
FX = FUNC ( X ) 

IF (KOUNT .GE. 60) WRITE (6,1004) KOUNT , X , FX , FOFX 
1004 FORMAT (1HL*9H KOUNT ,G16.9,9H X ,G16.9,9H FX 

19H FOFX , G 1 6 . 9 ) 

IF ( ABS ( FX-FOFX ) .LE. DELI) RETURN 
IF (KOUNT .EQ. 75) GO TO 1002 
IF ( A* ( FX-FOFX ) .LT. 0.) GO TO 1001 
XXO = X 

X = ( X + XX2 ) /2 . 

GO TO 1000 

1001 XX2 = X 

X = ( XXO+X ) / 2 . 

F2 = FX 
GO TO 1000 


, G 1 6 . 9 , 


1002 WRITE (6,1003) 

1003 FORMAT <//30X,62H75 
1ING TO A ROOT) 

RETURN 


ITERATIONS HAVE BEEN PERFORMED WITHOUT CONVERG 


1005 WRITE (6,1006) FOFX 

1006 FORMAT ( / / 1 OX , 7HF ( X ) 
RETURN 


,G16.9,31H IS OUTSIDE OF SPECIFIED LIMITS) 


1007 XI = XO 
RETURN 


1008 XI = X2 
RETURN 


END 
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$ I BFTC STARTT LIST 


SUBROUTINE START ( VL OW , F L 0 , VUP , FUP , V I MAX ) 

C OM MON /FAC TOR /PERM, SAME, GAM, GA MM 1,GAMP1, SSMLOW, SSMUP, REC ON V,BINTGR 
EXTERNAL CF AC T , QFAC T , FR AT , FOFRS , FKMAX 

XO = 1 ./PERM 

X2 = 0.999999999 

X I NTL = ( XO + X2 ) / 2 . 

IF ( VLOW .EO. 0.0) GO TO 70 

CALL ROOT (X0,X2, XINTL, FLU, FOFRS, RLOW) 

GO TO 71 

70 RLOW = 1.0 

71 SSMLOW = l./RLOW 

CALL ROOT < X0,X2 , XINTL , FLIP, FOFRS, RUP) 

SSMUP = 1 . /RUP 

IF (SSMLOW .EO. SSMUP) GO TO 40 

FKMAX(X) IS LINEAR IN A NEIGHBORHOOD OF X WHEN X IS 5UCH THAT FKMAX(X)=0 
USE GOOD INITIAL ESTIMATE PLUS LINEARITY TO FIND X SUCH THAT FKMAX(X)=0 

XINTL = ( l./PERM)*SORT ( SSMLOW/SSMUP ) 

XO = XINTL - 0.005 
FO = FKMAX ( XO , 0 ) 

X2 = XINTL + 0.001 
F2 = FKMAX ( X2 , 0 ) 

SLOPE = ( F2 - FO ) / ( X2 - XO) 

FINTL = FKMAX ( X INTL ,0 ) 

YINCPT = FINTL - SLOPE* X I NT L 
X A MK = -YINCPT/SLOPE 
FOF X = FKMAX ( XAMK , 1 ) 

IF (ABS(FOFX) .GT. 0.00009) WRITE (6,60) FOFX , XAMK 
60 FORMAT ( / / 2 9 X , 3 5H S E AR CH FOR ROOT FAILED F(X) = ,G16.9,7H X 

1= , G1 6 . 9 ) 

SAME = ( XAMK/SSMLOW) *( XAMK/SSMLUW ) 

C = 1. - PERM# ( GAMP1 ** ( 1 . / ( GAM— 1 . ) ) ) $ ( SSMUP/ ( SSMUP- SSMLOW ) ) *X AMK# 
1BINTGR 

CINTGR = SIMPSKSSMLOW, SSMUP, OFACT.K) 

0 = ( SSMLOW*SSMUP/ ( SSMUP-SSMLOW ) ) *C INTGR 
RATIO = 0/(1. - C) 

GO TO 50 

40 XAMK = l./PERM 

RATIO = SSMUP#SSMUP#OFACT ( SSMUP ) 

C = 0.0 
0 = 0.0 
50 XO = 1.0 
X2 = PERM 

X I NTL = ( XO + X 2 ) / 2 . 

CALL ROOT (X0,X2, XINTL, RATIO, FRAT,SSMIAX) 

VIMAX = ( 3 . 141 59265/4. )*( PERM-1 . ) + ( PERM/2 . ) *ARS I N ( 2 . #GAMM1 # 

1 SSM1 AX*SSMI AX - GAM) + 0 . 5# ARS I N ( 2 . *GAMP 1 / ( SSM I AX#SSM I A X ) - GAM) 
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VIMAX = VIMAX*RECONV 
VLOW = VLOW*R ECONV 
VUP = VUP^RECONV 

WSTAR = ( ( 1 . /GAMP1 ) =** ( GAMP1 / ( 2 .*GAMM1 ) ) )*CI NTGR 


10 


90 


20 


FDR SUPERSONIC STARTING) 


WRITE (6,10) 

FORMAT ( / /48X , 36HC ALCUL AT I ONS 
WRITE (6,90) WSTAR 

FORMAT ( /50X,24HWE IGHT-FLOW PARAMETER = ,F9,9) 

WRITE (6,20) XAMK,C,0,SSMIAX co , Ky 

FORMAT ( /20X , 10HK* ( MAX ) = ,F9.4,6X,5H C - ,F9.4,5X,5H 0 - ,F9. , , 

113HM*( I (MAX) ) = , F9 » 4 ) 

WRITE (6,30) VIMAX, VLOW, VUP, GAM 
30 FORMAT (/4X»38HTHE MAXIMUM DESIGN VALUE FOR V(IN) IS 

1 WHEN V(LOW) IS , F9 .4,1 6H DEG, V(UP) IS ,F9.4,16H 

2 ,F7.4) 

RETURN 
END 


,F9.4,21H DEG 
DEG, GAMMA IS 


$ I BFTC MSS LIST 

SUBROUTINE MSSTAR (M,N,VSSTAR) 


COMMON/ FAC TOR /PERM , SAME, GAM, GAMM1,GAMP1 , SSMLOW , SSMUP , RECON V , D 
REAL M , MS 
EXTERNAL ADSTR 


A = 0.785398162*(PERM -1.) 

B = 0.5*PERM 
C = GAM -1. 

0 = GAM + 1. 

I F ( N .NE. 0) GO TO 1 
MS = ADSTR ( M ) 

IF (MS .LT. 1.) GO TO 3 
GO TO 2 

1 XO = 1. 

X2 = PERM 

XINTL = ( XO + X2 ) /2. 

FOF X = M 

CALL ROOT (X0,X2, XINTL, FOFX, ADSTR, MS 

2 SQRDMS = MS*MS 

VSSTAR = ( A + B*AR SI N ( C* SQRDMS -GAM ) 
1 RECONV 
RETURN 


+ 0.5*ARSIN( D/SORDMS-GAM ) )4 


3 VSSTAR = 0. 
RETURN 


END 
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$ I BFTC SIMPS LIST 

FUNCTION SIMPS1 ( XM IN , XMAX , FUNC 1 , K ER ) 

DIMENSION V ( 200 ) * H < 200 ) ,A(200) , B ( 200 ) , C ( 200 ) , P ( 200 ) , F ( 200 ) , NE ( 2 00 ) 
EOOIVALENCE ( E , NE ) , ( T E ST , N T E ST ) 

T=3.0E-5 
V( I ) = XM I N 

H( 1 > =0 . 5* ( XMAX-XMIN) 

A ( 1 )=FUNC1 ( XMIN) 

B ( 1 ) =FUNC I ( XMIN+HI 1 ) ) 

C(l) =FUNC 1 ( XMAX ) 

P( 1 )=H( 1 )*(A( l)+4.0*d( 1 )+C< 1 ) ) 

E ( 1 ) =P ( 1 ) 

ANS=P ( 1 ) 

N = 1 

FRAC=2.0*T 

1 FR AC =0 . 5*FR AC 

2 TEST=ABS(FRAC#ANS) 

K = N 

3 DO 7 1=1, K 

A IF (NTEST-IABS(NE( I) ) ) 5,5,7 

5 N = N+l 

V< N » = V ( I ) + HI I ) 

H(N)=0.5*H< I ) 

A ( N > =B ( I ) 

B ( N ) =F UNC I ( V(N)+H(N) ) 

C ( N ) =C ( I ) 

P(N)=H(N)*(A(N)+4.0*B(N)+C(N) ) 

0=P( I ) 

H ( I ) =H ( N ) 

B ( I ) =FUNC 1 IVin+HlIH 
C ( I ) =A ( N ) 

P( I )=H( I ) * ( A ( I )+4.0*B ( I )+C( I ) ) 

Q=P( I ) +P ( N ) — 0 
ANS=ANS+0 
E ( I ) =0 
E ( N ) =0 

6 IF ( N-200 ) 7,13,13 

7 CONTINUE 

8 IF (N-K) 9,9,2 
90= 0.0 

10 DO 11 1=1 ,N 

1 1 0=Q+E ( I ) 

12 IF (ABS(Q)-T*ABS(ANS) ) 14,14,1 

13 K ER=K ER + 1 

14 ANS=0 . 0 

15 DO 16 I = 1 , N 

16 ANS=ANS+P ( I ) 

SIMPS1=( ANS+O/ 30.0)/ 3.0 

17 RETURN 

END 
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tIBFTC ALLFUN 

FUNCTION ALFUNC (A,R,Y) 


COMMON/ EXPAL F/GA ME 
COMMON/ FAC TOR /PERM 


, SAME , GAM, GAMM1 , GAMP 1 , SSMLOW , SSMUP , 


RECONV »GRTY 


ALFUNC = ( 1 • / Y ) * ( ( A - B*Y*Y ) **GAME XP ) 


RETURN 

END 


$ 1 BFTC BAKE 

FUNCTION CFACT (Y ) 

COMMON/FACTOR/PERM, SAME, GAM, GAMM1 , GAMP 1 , SSMLOW, SSMUP, RECONV, GRTY 
EXTERNAL ALFUNC 

CFACT = ALFUNC ( 1 .» SAME , Y ) 

RETURN 

END 


$ I BFTC CHARL 

FUNCTION OFACT (Y) 

COMMON/ FACTOR/ PERMt SAME t GAM, GAMM1 , GAMP 1 » SSMLOW , SSMUP t RECONV » GRT Y 

EXTERNAL ALFUNC 

OFACT = ALFUNC ( GAM PI f GAMM1 , Y ) 

RETURN 

END 


$ I BFTC DOGG 

FUNCTION FRAT(Y) 

COMMON/ F AC TOR/ PERM » SAME t GAM tGAMMl ,GAMP 1 , SSMLOV* 

EXTERNAL OFACT, ALFUNC 

F rAT = ( Y**( GAM/GAMM1 )) *OF ACT ( Y )/ ALFUNC ( ~GAMM1 


, SSMUP, RECONV, GRTY 
,-GAMPl ,Y ) 


RETURN 

END 


SIBFTC FELI 

FUNCTION FOFRS (X) 


COMMON/FACTOR/PERM 
double precision X 


SAME, GAM, GAMM1, GAMP l,SSMLOW, SSMUP, RECONV,GRTY 


ARG1 = 2.*GAMM1/(X*X) - GAM 
ARG2 = 2.*GAMP1*X*X - GAM 

, I ^ BSIW1) * GT * 1,0 *° R ‘ ab S«A«G2) .GT. 1.0) WRITE 

I f A K l? c 

1 FORMAT I//14X.61HARGUMENT OF ARCSIN IS OUTSIDE DOMAIN OF 
IN ARG1 = ,G16.9 ; 1 H ARG2 = ,G16.9) 


(6,1) ARG 1 
DEF INITIO 


FOFRS = PERM*ARS I N ( ARG1 ) + ARS I N ( ARG2 ) 


RETURN 

END 


SIBFTC GERT 

FUNCTION FKMAX ( Y , L ) 


C OMMON/ FAC TOR/PERM , SAME, GAM, GAMM1, GAMP 1 
EXTERNAL AL FUNC ,CFACT 


SSMLOW,SSMUP,RECONV,BINTGR 


SAME = (Y/SSMLOW)*(Y/SSMLOW) 

K = 0 

FKMAX = S I MPS 1 ( SSMLOW , SSMUP , CF ACT,K ) 

IF (K .EO. 1 ) WRITE (6,1 ) 

1 FORMAT ( / 1 OX , 2 6HFA I LURE TO INTEGRATE CFACT) 

IF (L .EO. 1) B I N TGR = FKMAX 

FKMAX = FKMAX + SSMUP*CFACT < SSMUP) - ALFUNC ( 1 . , Y*Y, 1 . ) 

RETURN 

END 
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$ I BF TC Sl'ARM 

FUNCT 1 1 IM ADSTR ( MS TAR ) 

COMMON/ F ACTOR/ PERM, SAME, GAM, GAMM1 , GAMP 1 , SSMl UW , S SMUP , R ECONV , 0 
KFAL MSTAR,M 


C = GAM - 1. 

F = C/O 
G = C/GAM 
H = GAM/ 0 
M = MST AR*MST AR 

ADS T R = PERM*SORT( ( 1 • - ( 1 • _ E #M ) * ( ( 1 . +0 . 5 * 


H*M ) / ( 1 .-F*M ) ) ) **G ) ) ) 


RETURN 

END 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 6, 1967, 
128-31-02-25-22. 
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